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Introduction: 


Enhanced  MRI  imaging  techniques  called  restriction  spectrum  imaging  (RSI)  has  been  utilized  to 
identify  glioblastoma  multiforme  (GBM)  brain  tumors  (1).  The  output  from  these  techniques  has 
been  termed  the  tumor  ‘RSI  cellularity  index’  and  provides  significantly  greater  accuracy  in 
distinguishing  brain  tumor  from  normal  tissue  than  traditional  imaging  measures.  The  RSI  tumor 
signal  was  noted  to  be  10-fold  greater  than  conventional  imaging  measures  (1-2).  In  addition, 
these  imaging  techniques  are  corrected  for  spatial  distortions  that  allows  for  more  accurate 
localization  of  the  tumor  signal  (3).  In  prostate  cancer,  increasing  Gleason  score  correlates  with 
loss  of  normal  gland  formation,  loss  of  peripheral  gland  tubular  structure  and  increased 
cellularity.  Therefore,  we  hypothesize  that  the  RSI  cellularity  index  will  correlate  with  higher 
tumor  grade  based  on  Gleason  score  and  will  provide  significantly  greater  accuracy  in 
discriminating  aggressive  tumor  from  benign  and  indolent  lesions  when  compared  to  current 
imaging  techniques.  Based  on  its  prior  ability  to  discriminate  between  brain  tumors  and  adjacent 
normal  tissue  with  high  accuracy  and  high  signal  to  noise,  we  hypothesize  that  the  RSI  cellularity 
index  will  perform  similarly  in  prostate  tissue.  We  propose  a  non-invasive  (without  endorectal 
coil  or  IV  contrast),  rapid  imaging  test  based  on  these  procedures.  The  goal  of  the  current 
research  is  to  test  whether  that  objective  is  feasible. 

Keywords: 

DCE  =  dynamic  contrast  enhancement 
DWI  =  diffusion  weighted  imaging 
ADC  =  apparent  diffusion  coefficient 
RSI  =  restriction  spectrum  imaging 
Cl  =  cellularity  index 
EPI  =  echo  planar  imaging 
Prostate  cancer  screening 
Prostate  cancer  staging 
Non-invasive  imaging 
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Overall  Project  Summary: 

Specific  Aim  1 : 

All  aspects  of  Specific  Aim  1  have  been  accomplished.  Existing  pulse  sequences  at  UCSD  have 
been  modified  for  signal  optimization  in  patients.  Post-processing  techniques  for  displaying  and 
analyzing  RSI  data  have  been  developed  and  optimized.  The  whole  mount  sectioning  protocol 
for  prostate  specimens  at  UCSD  has  been  optimized.  Local  IRB  approval  was  obtained  under 
HRPP#071983,  which  was  also  accepted  by  DoD. 

Specific  Aim  2: 

Major  Tasks  2  and  3  include  data  collection  and  analysis  and  is  ongoing.  MR  imaging  with  an 
endorectal  coil  for  85  patients,  without  an  endorectal  coil  for  10  patients,  and  TIFFs  of  whole- 
mount  pathology  slides  for  28  patients,  have  been  anonymously  transferred  to  UCSD  from 
UCLA  for  analysis.  In  parallel  we  have  been  collecting  imaging  and  pathology  data  at  UCSD 
without  an  endorectal  coil  using  similar  imaging  and  pathology  techniques.  All  imaging 
evaluated  from  UCSD  and  UCLA  have  RSI  techniques  incorporated.  As  noted  below  we  have 
imaged  182  patients  at  UCSD  with  17  whole  mount  pathology  sections.  We  have  currently 
collected  roughly  33%  of  our  goal  of  150  patients  with  whole  mount  pathology  sections. 
Analysis  of  these  data  has  so  far  resulted  in  one  published  manuscript  (White  et  al,  Cancer 
Research  2014),  one  accepted  manuscript  (Rakow-Penner,  et  al.)  and  one  submitted  manuscript 
(Liss  et  al.).  In  summary,  all  tasks  and  subtasks  of  Aim  2  have  been  partially  perfonned. 
Completion  will  follow  after  all  patients  have  been  imaged  over  the  next  few  months. 


Specific  Aim  3: 

Major  Task  4  and  5  includes  data  collection  and  analysis  and  is  ongoing.  Imaging  using  RSI  has 
been  implemented  as  part  of  the  patient  standard  of  care  at  UCSD,  both  for  screening  and  for 
pre-operative  MRI.  All  patients  undergoing  pre-operative  MRI  at  UCSD  are  scanned  using  the 
modified  imaging  protocol  without  an  endorectal  coil,  resulting  in  182  scans  at  UCSD  during  the 
covered  grant  period.  Prostate  tissue  is  obtained  from  each  prostatectomy  patient  and  processed 
using  the  optimized  whole  mount  sectioning  protocol.  Seventeen  patients  at  UCSD  underwent 
prostatectomy  with  subsequent  whole-mount  pathology  during  the  covered  grant  period.  We  are 
continuing  to  collect  imaging  and  pathology  data.  Our  stated  goal  in  the  SOW  was  to  collect 
imaging  and  whole  mount  path  data  from  75  patients  by  month  17.  We  expect  to  reach  this  goal 
as  patients  with  imaging  over  the  last  year  will  be  finally  going  to  prostatectomy.  We  can 
expect  to  perform  a  similar  analysis  as  in  Aim  2  in  this  subset  of  patients  that  were  imaged 
without  an  endorectal  coil.  In  summary,  all  tasks  and  subtasks  of  Aim  3  have  been  partially 
performed.  Completion  will  follow  after  all  patients  have  been  imaged  over  the  next  few 
months. 
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Key  Research  Accomplishments: 

1)  We  have  identified  spatial  distortion  in  conventional  diffusion  weighted  imaging  and,  using 
RSI,  we  are  able  to  correct  for  it.  This  table  summarizes  the  results  for  28  patients.  The  results 
show  that  distortions  of  tumor  signal  can  be  as  great  as  1.3  cm  if  left  uncorrected. 


Table  2.  Distortion  distance  of  the  whole  prostate  and  tumor  regions  of  interest  due  to  B0 
inhomogeneity 

Mean  (standard  deviation) 
(mm) 

Maximum  distortion 
(mm) 

Whole  prostate  distortion 

3.2  (2.2) 

12 

Tumor  ROI  distortion 

3.2  (2.4) 

13 

Rakow-Penner  et  al,  in  preparation 


2)  As  a  result  of  the  correction  for  distortion,  RSI  has  improved  our  detection  of  extra  prostatic 
extension  by  more  accurately  localizing  tumor  location. 


MR  imaging  and  whole  mount  pathology  for  7  prostate  cancer  patients.  Column  1:  T2  -  weighted  images,  column 
2:  Perfusion  K^s  maps,  column  3:  standard  apparent  diffusion  coefficient  (ADC)  maps  (b  =  0,  100,  400,  and  800 
s/mm2),  column  4:  RSI  z-score  maps,  column  5:  whole  mount  pathology  with  tumor  and  area  of  EPE  identified. 
Color  bar  represents  zscores  from  0  to  7  for  the  RSI  maps.  (Rakow-Penner  et  al,  accepted) 
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3)  RSI  imaging  has  also  been  implemented  as  part  of  the  work  flow  for  defining  tumor  regions  of 
interest  for  targeted  biopsy  using  Dynacad  and  UroNav  MR-fused  ultrasound  targeted  biopsy 
software.  Due  to  increased  conspicuity  of  the  RSI  signal  compared  to  perfusion  (DCE)  and 
conventional  DWI,  as  well  as  correction  of  spatial  distortions,  we  are  able  to  localize  small 
tumors  with  accuracy  as  illustrated  below.  While  this  was  not  an  Aim  of  the  DoD  proposal,  this 
is  certainly  a  noteworthy  achievement. 


Perfusion 

T2  Perfusion  Time  Series  RSI  MIP 


Images  from  three  patients  who  underwent  targeted  biopsy.  The  first  column  demonstrates  the  T2  map  with  the  ROI 
drawn  to  guide  the  MR-fused  ultrasound  guided  biopsy.  The  second  column  displays  the  Dynacad  perfusion  maps, 
where  red  signifies  rapid  intense  enhancement  with  washout,  green  signifies  rapid  intense  enhancement  with  plateau 
washout  and  blue  signifies  gradual  enhancement.  The  third  column  demonstrates  the  perfusion  curves  for  the 
location  identified  by  the  crosshair  in  the  other  images.  The  fourth  column  illustrates  a  heatmap  based  on  the  RSI 
overlaid  on  T2  images,  with  red  regions  correlating  to  the  areas  of  greatest  restricted  diffusion.  The  fifth  column 
shows  maximum  intensity  projections  of  the  prostate  boundary  (green)  and  ROI  for  targeting  (yellow).  A)  Patient 
with  three  negative  previous  biopsies.  The  targeted  biopsy  found  Gleason  4+3  cancer.  B)  Patient  with  previous 
incomplete  radiation  therapy  for  response  to  low-grade  cancer.  Targeted  biopsy  found  Gleason  4+4  cancer.  C) 
Patient  with  previous  biopsy  demonstrating  low-grade  cancer.  Targeted  biopsy  found  Gleason  3+4  cancer. 
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4)  Quantitatively,  we  have  shown  RSI  to  be  at  least  as  good  as  ADC  in  detecting  high- 
grade  cancer.  Our  results  suggest  that  the  RSI  z-score  may  accurately  predict  Gleason 
grade.  The  RSI  z-score  is  a  quantitative  measure  of  cellularity. 


T2  ADC  RSI  Histopath 


Representative  images  showing  RSI  z-score  maps  across  Gleason  scores.  The  y-axis  shows  the  pathologic  Gleason 
score  with  the  x-axis  designating  the  MRI  sequence.  The  last  column  displays  the  whole  mount  pathology  with  the 
corresponding  cancer  region  of  interest  circled  in  black.  The  star  3+3  pathologic  figure  represents  the  lesion  with 
3+3  prostate  cancer  while  the  other  lesions  are  4+3.  White  arrows  in  the  4+3  and  5+4  RSI  indicate  regions  where  the 
ADC  is  dark  and  would  be  suspicious  for  tumor,  but  the  RSI  reveals  these  are  regions  with  signal  void  and  are 
unlikely  to  be  cancerous.  Color  bar  represents  z-scores  from  -2.8  to  7  for  the  RSI  maps.  (Liss  et  al,  submitted) 


RSI  Z  Score  Per  Region  of  Interest 


RSI  z-score  value  grouped  by  pathologic  Gleason  score  for  28  patients:  The  y  -axis  represents  the  RSI  z-score 
derived  from  a  given  region  of  interest.  The  x-axis  demonstrates  the  pathological  Gleason  scores  in  increasing  levels 
of  aggressiveness  from  left  to  right.  Each  data  point  represents  one  region  of  interest  corresponding  to  a  location  on 
the  whole  mount  prostatectomy  specimen  contoured  by  a  GU  pathologist  (Liss  et  al,  submitted). 
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Conclusion: 

•  We  have  identified  the  distortions  in  RSI  and  conventional  diffusion  weighted  imaging 
and  corrected  for  them  in  our  postprocessing  stream.  This  has  increased  our  ability  to 
detect  extraprostatic  extension  and  helped  to  increase  the  accuracy  of  targeted  biopsies. 

•  We  have  optimized  and  implemented  a  non-invasive  imaging  measure  of  Gleason  grade 
(RSI  z  -score)  that  is  at  least  as  good  as  standard  ADC. 

•  All  RSI  cases  scanned  at  UCSD  have  been  performed  without  an  endorectal  coil  and  with 
sequences  of  less  than  10  minutes,  suggesting  the  feasibility  of  using  this  protocol 
clinically  on  a  non-invasive  screening  basis. 
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and  Applications  of  Restriction  Spectrum  Imaging 
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Natalie  M.  Schenker-Ahmed1,  Hauke  Bartsch1,  Rebecca  Rakow-Penner1,  Dominic  Holland1,  Ahmed  Shabaik3, 
Atle  Bjornerud4,  Tuva  Hope5,  Jona  Hattangadi-Gluth6,  Michael  Liss7,  J.  Kellogg  Parsons7,  Clark  C.  Chen8, 
Steve  Raman9,  Daniel  Margolis9,  Robert  E.  Reiter10,  Leonard  Marks10,  Santosh  Kesari11,  Arno  J.  Mundt6, 
Christopher  J.  Kane7,  Bob  S.  Carter8,  William  G.  Bradley1,  and  Anders  M.  Dale1,11 


Abstract 

Diffusion-weighted  imaging  (DWI)  has  been  at  the  forefront  of  cancer  imaging  since  the  early  2000s.  Before  its 
application  in  clinical  oncology,  this  powerful  technique  had  already  achieved  widespread  recognition  due  to  its 
utility  in  the  diagnosis  of  cerebral  infarction.  Following  this  initial  success,  the  ability  of  DWI  to  detect  inherent 
tissue  contrast  began  to  be  exploited  in  the  field  of  oncology.  Although  the  initial  oncologic  applications  for  tumor 
detection  and  characterization,  assessing  treatment  response,  and  predicting  survival  were  primarily  in  the  field 
of  neurooncology,  the  scope  of  DWI  has  since  broadened  to  include  oncologic  imaging  of  the  prostate  gland, 
breast,  and  liver.  Despite  its  growing  success  and  application,  misconceptions  about  the  underlying  physical  basis 
of  the  DWI  signal  exist  among  researchers  and  clinicians  alike.  In  this  review,  we  provide  a  detailed  explanation  of 
the  biophysical  basis  of  diffusion  contrast,  emphasizing  the  difference  between  hindered  and  restricted  diffusion, 
and  elucidating  how  diffusion  parameters  in  tissue  are  derived  from  the  measurements  via  the  diffusion  model. 
We  describe  one  advanced  DWI  modeling  technique,  called  restriction  spectrum  imaging  (RSI).  This  technique 
offers  a  more  direct  in  vivo  measure  of  tumor  cells,  due  to  its  ability  to  distinguish  separable  pools  of  water  within 
tissue  based  on  their  intrinsic  diffusion  characteristics.  Using  RSI  as  an  example,  we  then  highlight  the  ability  of 
advanced  DWI  techniques  to  address  key  clinical  challenges  in  neurooncology,  including  improved  tumor 
conspicuity,  distinguishing  actual  response  to  therapy  from  pseudoresponse,  and  delineation  of  white  matter 
tracts  in  regions  of  peritumoral  edema.  We  also  discuss  how  RSI,  combined  with  new  methods  for  correction  of 
spatial  distortions  inherent  in  diffusion  MRI  scans,  may  enable  more  precise  spatial  targeting  of  lesions,  with 
implications  for  radiation  oncology  and  surgical  planning. 

See  all  articles  in  this  Cancer  Research  section,  "Physics  in  Cancer  Research." 
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Introduction 

Diffusion-weighted  imaging  (DWI)  is  increasingly  used  as  an 
imaging  biomarker  for  the  detection  and  characterization  of 
primary  and  metastatic  brain  tumors  as  well  as  for  prognosti¬ 
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cation  and  monitoring  treatment  response  in  this  patient  pop¬ 
ulation.  The  unique  ability  of  this  technique  to  probe  the 
underlying  structure  of  brain  tissue  at  a  cellular  level  makes  it 
well  poised  to  answer  questions  about  tumor  biology,  as  well  as 
the  microstructure  of  peritumoral  white  matter.  The  greatest 
enthusiasm  to  date  has  surrounded  the  ability  of  DWI  to 
estimate  tumor  cellularity  on  the  basis  of  quantitative  images 
of  diffusion,  namely  apparent  diffusion  coefficient  (ADC)  maps. 
In  particular,  numerous  investigators  have  shown  a  negative 
correlation  between  tumor  ADC  values  and  cellularity  (1-3), 
which  is  commonly  attributed  to  increased  restricted  diffusion 
imposed  by  tumor  cells.  Following  this  logic,  ADC  has  been  used 
with  varying  success  to  characterize  and  grade  primary  and 
metastatic  brain  tumors  (3-9),  to  assess  tumor  response  to 
therapy  (10-12),  and  to  predict  survival  in  patients  with  malig¬ 
nant  tumors  (13-17).  However,  confusion  about  the  true  bio¬ 
physical  basis  of  the  DWI  signal  abounds,  and  awareness  of  the 
limitations  of  the  existing  approaches  has  led  to  the  develop¬ 
ment  of  advanced  DWI  methods  that  extend  beyond  the  tensor 
model.  These  methods  include  high-angular  diffusion  imaging 
techniques,  such  as  diffusion  spectrum  imaging  (DSI;  ref.  18)  and 
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Q-ball  imaging  (19),  as  well  as  methods  to  probe  non-Gaussian 
diffusion,  including  biexponential  (20, 21),  stretched  exponential 
(22),  and  kurtosis  imaging  (23),  and  methods  to  study  perfusion- 
related  effects  at  low  ft-values,  such  as  intravoxel  incoherent 
motion  (IVIM;  ref.  24).  Although  these  techniques  often  provide  a 
better  characterization  of  tissue  architecture  than  traditional 
models,  the  relationship  of  these  measures  to  the  underlying 
pathophysiology  of  tumors  is  largely  unknown. 

Advances  in  magnetic  resonance  (MR)  technology  are  now 
creating  even  more  possibilities,  increasing  excitement  in  the 
field  by  broadening  the  potential  applications  of  DWI  within 
cancer  diagnosis  and  treatment.  Newer  scanners  with  improved 
gradient  performance  allow  for  higher  ft- values  with  shorter 
diffusion  and  echo  times.  In  addition,  improved  coil  technology 
and  higher  field  strength  provide  better  signal-to-noise  ratio 
(SNR).  With  these  MR  advancements,  it  has  become  possible  to 
develop  new,  innovative  DWI  methods  that  provide  more  direct 
measures  of  tumor  cellularity  by  leveraging  the  intrinsic  con¬ 
trast  of  tumor  cells  relative  to  other  tissues.  Such  direct  mea¬ 
sures  address  a  top  challenge  posed  by  the  NCI  to  develop  in  vivo 
imaging  methods  that  portray  tumor  "cytotypes,"  i.e.,  imaging 
methods  that  can  probe  the  identity,  quantity,  and  location  of 
different  cells  that  make  up  a  tumor  and  its  microenvironment 
(NIH  RFA-CA- 13-020).  However,  the  successful  application  of 
such  methods  requires  not  only  advanced  MR  technology  and 
an  appreciation  of  the  clinical  challenges  in  neurooncology, 
but  also  an  exquisite  understanding  of  the  physical  basis  of  the 
DWI  signal  and  its  current  limitations. 

The  purpose  of  this  review  is  to  provide  a  detailed  explanation 
of  the  biophysical  basis  of  diffusion  contrast  and  to  demon¬ 
strate  what  is  known  about  how  it  reflects  tissue  microstructure 
in  the  context  of  key  clinical  dilemmas  in  neurooncology.  We 
begin  with  a  basic  introduction  of  how  diffusion  measurements 
are  derived  from  DWI  and  then  follow  with  a  review  of  the 
biophysical  properties  of  water  diffusion  in  tissue.  We  then 
describe  how  diffusion  parameters  in  tissue  are  derived  from  the 
measurements  via  fitting  of  the  diffusion  model  to  observed 
data.  For  this  review,  we  focus  on  one  advanced  DWI  modeling 
technique  called  restriction  spectrum  imaging  (RSI),  which 
provides  a  general  framework  for  estimating  tissue  properties 
from  DWI  data  and  addresses  the  NCI  challenge  for  more  direct, 
in  vivo  imaging  of  tumor  cells.  We  also  address  the  power  of 
multispectral  imaging,  highlighting  the  importance  of  combin¬ 
ing  advanced  DWI  methods  such  as  RSI  with  other  imaging 
modalities  (i.e„  MRI  perfusion,  PET)  to  optimize  the  detection 
and  monitoring  of  brain  tumors.  Although  the  primary  focus  of 
this  review  is  on  the  application  of  RSI  to  neurooncology,  the 
same  methods  are  equally  applicable  to  other  applications  of 
oncologic  imaging,  such  as  prostate,  liver,  and  breast  cancer, 
and  we  briefly  present  some  preliminary  data  in  prostate  cancer 
at  the  end  of  "clinical  applications." 

Principles  of  Diffusion  MRI 

The  diffusion  experiment 

While  early  diffusion  measurements  using  nuclear  magnetic 
resonance  (NMR)  in  biologic  tissue  were  made  in  the  1960s  and 
1970s,  it  was  not  until  the  mid  1980s  that  the  basic  principles  of 
DWI  emerged  (25-27).  All  MRI  techniques,  including  DWI,  rely 


on  the  fundamental  principle  that  water  hydrogen  nuclei 
become  magnetized  when  placed  in  a  strong  static  magnetic 
field.  MRI  contrast  is  formed  by  perturbing  this  magnetization 
using  electromagnetic  waves  and  then  allowing  the  magneti¬ 
zation  to  "relax"  back  to  an  equilibrium  state.  Tissue  relaxation 
times — namely,  the  spin-lattice  ( 7', )  and  spin-spin  ( T2)  relax¬ 
ation  times — form  the  fundamental  basis  of  soft  tissue  contrast 
and  anatomic  imaging  with  MRI.  DWI  is  an  MRI  technique  that 
adds  additional  sensitivity  to  the  intrinsic  random  thermal 
displacements  of  water  molecules  that  take  place  during  the 
MRI  experiment  as  part  of  the  natural  diffusion  process.  A  pair 
of  pulsed  magnetic  field  gradients  are  turned  on  and  off  in 
succession  (with  duration  8)  to  magnetically  "label"  and  "refo¬ 
cus”  the  spin  phase  of  hydrogen  nuclei,  depending  on  where 
they  reside  physically  in  the  gradient  field  (Fig.  1).  The  second 
refocusing  pulse  is  applied  some  finite  time  A  after  the  first 
pulse  in  an  effort  to  realign  the  spin  phases  of  stationary  nuclei. 
Thus,  any  residual  spin  phase  left  after  the  application  of  the 
refocusing  pulse  can  be  attributed  to  the  diffusion  of  water 
along  the  orientation  of  the  field  gradient  during  the  effective 
diffusion  time  of  the  experiment  TA,  defined  here  as  A  —  (8/3). 
The  net  phase  dispersion  due  to  diffusion  causes  an  attenu¬ 
ation  of  the  measured  signal  and  a  decrease  in  the  voxel 
intensity.  Stejkal  and  Tanner  (28)  were  some  of  the  pioneers 
of  the  pulsed  field  gradient  experiment  and  provided  the 
mathematic  framework  to  relate  the  diffusion  coefficient  to 
the  experimental  variables.  To  date,  the  Stejkal  Tanner  pulse 
sequence  remains  the  gold  standard  method  for  measuring 
diffusion  in  the  clinic  and  forms  the  basis  of  quantitative 
mapping  of  tumor  cellularity  with  ADC,  as  discussed  below. 

Biophysics  of  water  diffusion  in  tissue 

In  this  section,  we  provide  a  basic  review  of  the  three 
principal  physical  modes  of  diffusion  in  tissue:  free,  hindered, 
and  restricted. 

Free  diffusion.  Free  water  diffusion  describes  the  random 
(Brownian)  motion  of  water  molecules  due  to  thermal  agita¬ 
tion,  in  the  absence  of  any  obstacles.  The  displacement  dis¬ 
tribution  of  free  water  molecules  is  time  dependent  and 
Gaussian  and  obeys  a  statistical  law  established  by  Einstein 
in  1905.  Along  a  single  direction  in  space,  the  average  molecular 
excursion  of  water  molecules  is  proportional  to  the  square  root 
of  diffusion  time  TA.  Mathematically,  this  can  be  expressed  as 
s  =  (2 DTd)1/2,  where  s  is  the  root-mean-squared  distance  and  D 
is  the  diffusion  coefficient  (Fig.  2A).  For  free  water  at  brain 
temperatures  (37°  C),  the  diffusion  coefficient  is  approximately 
3  |lm2/ms  (29),  which  translates  to  a  distance  of  approximately 
17  |lm  in  50  milliseconds.  In  brain  tissue,  however,  water 
molecules  are  constantly  bouncing  off  and  interacting  with 
various  tissue  elements  such  as  cell  membranes  and  macro¬ 
molecules.  The  net  result  is  a  decrease  in  diffusion  mobility  of 
water  and  a  displacement  distribution  that  no  longer  follows  a 
single  Gaussian  distribution.  Generally  speaking,  tissue  ele¬ 
ments  impede  water  diffusion  both  through  diffusion  hin¬ 
drance  and  restriction.  Hindered  and  restricted  diffusion  are 
two  distinct  processes  that  result  from  fundamentally  different 
behavior  of  spins  within  the  intra-  and  extracellular  tissue 
compartments,  as  described  next. 
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Figure  1.  The  diffusion  experiment. 
Sensitivity  to  the  random  molecular 
displacements  (Brownian  motion) 
of  water  molecules  is  achieved 
through  the  use  of  two  magnetic 
field  gradient  pulses  with  amplitude 
G,  duration  <5,  and  separation  A. 
During  the  first  pulse,  the  initial 
positions  of  water  molecules 
(spins)  are  encoded  with  a  phase 
offset,  depending  on  their  spatial 
location  in  the  gradient  field.  The 
second  pulse  is  then  applied  after 
some  finite  delay  A  to  realign  the 
spin  phases.  In  this  way,  if  water 
molecules  diffuse  to  a  different 
physical  location  along  the  gradient 
field  direction,  refocusing  will  be 
imperfect  and  a  net  phase 
dispersion  will  result.  This  phase 
dispersion  causes  an  attenuation  of 
the  magnitude  signal  and  a 
decrease  (darkening)  of  the 
measurement  voxel  in  the 
reconstructed  image. 


Hindered  diffusion.  Hindered  diffusion  refers  to  the  delay 
of  passage  of  small  molecules  as  they  navigate  around  cellular 
obstacles,  as  in  brain  extracellular  space  (ECS;  Fig.  2B).  There  is 
a  long  history  (30)  of  quantifying  the  degree  of  hindrance  to 
diffusion  in  the  ECS  compared  with  that  of  a  free  medium  by  the 
tortuosity  X  (31),  defined  simply  as  the  square  root  of  the  ratio  of 
the  free  (or  intrinsic)  extracellular  diffusion  coefficient  Dt.MrA  to 
the  measured  or  ADC  for  extracellular  water  ADCextra:  X  = 
(Ajxua/ADCextra)172.  As  such,  tortuosity  simply  reflects  the 
degree  to  which  the  ECS  slows  diffusion  relative  to  free  water. 
Extensive  experimental  studies  using  real-time  iontophoresis 
(RTI)  using  the  cation  tetramethylammonium  (TMA)  indicate 
that  diffusion  through  tortuous  ECS  in  normal  brain  is  generally 
slowed  by  a  factor  of  about  2.6,  corresponding  to  X  of  about  1.6 
(31).  Although  there  may  be  several  factors  that  ultimately 
contribute  to  the  delay  of  passage  of  water  molecules  in  brain 
ECS  (such  as  the  composition  of  the  extracellular  matrix  and 
the  transient  trapping  or  binding  of  water  molecules  to  cell 
surfaces;  see  ref.  32),  perhaps  the  greatest  contribution  is  ECS 
geometry  (dictated  by  the  extracellular  volume  fraction  a  and 
shape  of  cells).  As  a  decreases,  due  to,  for  example,  cell  swelling 
or  greater  cell  packing  density  (cellularity),  water  molecules 
must  travel  more  circuitous  paths  around  cellular  obstructions 
and  tortuosity  increases.  The  maximum  theoretical  tortuosity 
due  to  geometry  can  be  quantified  as  the  ratio  of  the  distance 
though  the  center  of  a  spherical  cell  to  the  distance  around  the 
periphery  of  the  cell  or  Xg  =  n/2  =  1.57  (Fig.  2,  inset).  However, 
more  complex  simulations  over  a  wide  variety  of  packed  cellular 
objects  indicate  that  the  maximum  tortuosity  due  to  geometry 
is  no  greater  than  1.22  (33).  In  tumor  ECS  (or  "tumor  inter- 
stitium"),  both  tortuosity  and  volume  fraction  a  are  generally 
higher.  For  example,  in  anaplastic  astrocytomas  and  glioblas¬ 


tomas,  the  average  a  was  measured  to  be  as  high  as  0.47  to  0.49, 
with  X  around  1.67  to  1.77  (34).  Increased  a  has  been  associated 
with  necrosis  and  or  vasogeneic  edema,  whereas  increased  X 
may  result  from  either  astrogliosis,  commonly  observed  in 
tumor  tissue,  or  to  changes  in  the  extracellular  matrix 
(32,  34).  In  summary,  while  it  is  often  suggested  that  the  low 
ADC  observed  in  high-grade  tumors  results  from  increased 
packing  density  (cellularity)  of  cells  due  to  greater  hindrance 
imposed  on  extracellular  water  diffusion,  from  a  physical 
perspective,  cellular  crowding  can  only  mildly  reduce  the 
effective  diffusion  coefficient  in  tumors  (A.  increases  only  nar¬ 
rowly  with  decreasing  a).  Therefore,  it  stands  to  reason  that  a 
major  component  in  decreased  ADC  in  tumors,  especially  at 
high  Z>- values,  is  restricted  diffusion  within  the  cellular  com¬ 
partments  themselves,  described  next. 

Restricted  diffusion.  Restricted  diffusion  is  a  term  clas¬ 
sically  used  (35)  to  describe  the  trapping  of  water  molecules 
within  an  enclosed  compartment  (i.e„  as  defined  by  the  cell 
plasma  membrane)  such  that  the  net  distance  traveled  is  re¬ 
stricted  or  confined  by  the  compartment  dimensions. 
Although  the  term  "restricted  diffusion"  is  often  used  in  the 
clinical  literature  to  refer  to  any  reported  decrease  in  ADC,  the 
physics  of  restricted  diffusion  is  fundamentally  different  from 
hindered  or  free  water  diffusion.  Specifically,  for  restricted 
diffusion,  the  net  squared  displacement  of  water  molecules  is 
sublinear  in  time  and  therefore  non-Gaussian  (Fig.  2C).  In 
addition,  the  time  evolution  of  net  displacements  strongly 
depends  on  the  size  and  shape  of  the  restricting  compartment 
as  well  as  the  intrinsic  intracellular  diffusivity  Antra-  In  con¬ 
tradistinction  to  hindered  extracellular  water,  where  ADCextra 
is  independent  with  diffusion  time  Td,  the  effective  diffusion 
coefficient  for  intracellular  water  (or  ADCjntra)  decreases  with 
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Figure  2.  The  three  principal  modes  of  diffusion  in  tissue.  A,  free  diffusion;  in  free  water,  the  average  molecular  excursion  along  a  single  dimension  in  space  in 
terms  of  the  root-mean-squared  distance  s  increases  linearly  with  the  square  root  of  diffusion  time  s  =  (2DTd)1/2  with  a  slope  that  depends  on  the 
intrinsic  diffusivity  D.  B,  hindered  diffusion;  for  hindered  water  in  brain  ECS,  the  net  displacements  remain  linear  with  the  square  root  of  diffusion  time 
(i.e.,  Gaussian),  but  the  effective  diffusion  coefficient  D*  (or  ADC)  is  reduced  compared  with  D  due  to  tortuosity  of  the  ECS.  The  theoretical  maximum 
reduction  in  D*  (or  ADC)  that  can  be  expected  due  to  crowding  of  small  spherical  cells  in  the  ECS  is  given  by  the  tortuosity  limit  n/2  or  40%  (35).  C,  restricted 
diffusion;  in  restricted  intracellular  diffusion,  the  net  distance  traveled  by  water  molecules  is  limited  by  the  compartment  dimensions,  leading  to  a 
sublinear  time  evolution  of  the  net  squared  displacement  and  a  decreased  ADC.  The  ADC  of  restricted  intracellular  water  decreases  with  diffusion  time  as  a 
larger  proportion  of  the  spins  "bounce  off"  the  plasma  membrane. 


Td  over  a  finite  range  as  a  larger  proportion  of  the  spins  "bounce 
off  the  plasma  membrane. 

It  is  important  to  note  that  restricted  (or  non-Gaussian) 
diffusion  must  be  discussed  within  the  context  of  diffusion 
time,  permeability,  and  the  size  of  the  restricting  cellular 
compartments.  For  instance,  in  the  short-time  limit  (Td  <  1 
ms),  diffusion  is  largely  unrestricted  and  dictated  by  the 
intrinsic  diffusivity  of  the  medium  (ADCintra  ~  Antra),  except 
for  a  small  minority  of  spins  located  in  close  proximity  to  cell 
membranes  or  other  barriers.  Conversely,  in  the  long-time  limit 
(rd  >  1  s),  much  greater  than  the  average  residence  time  for 
water  within  intra-  and  extracellular  compartments,  diffusion 
is  dominated  by  exchange  and  can  be  accounted  for  by  a  single 
apparent  tensor,  reflecting  the  effective  medium  approxima¬ 
tion  of  the  tissue  (ADCintra  ~  ADCtissue)-  In  the  intermediate¬ 
time  regimen,  relevant  to  typical  clinical  DWI  acquisitions  (TA 
~  50-100  ms),  the  diffusion  time  is  short  relative  to  exchange 
between  intra-  and  extracellular  compartments  through  the 
plasma  membrane  (36),  but  long  enough  for  water  molecules  to 
repeatedly  come  in  contact  with  and  "bounce  off  the  plasma 
membrane.  In  this  case,  the  intracellular  spins  approach  a  "fill¬ 
up"  regimen,  where  the  displacements  of  spins  are  physically 
restricted  by  the  plasma  membrane  and,  therefore,  dictated  by 
the  size  and  shape  of  the  cellular  compartment.  Note  that  in  the 
intermediate-time  regimen,  diffusion  in  the  ECS  remains  hin¬ 


dered  and  behaves  according  to  a  classic  Gaussian  model, 
where  the  tortuosity  is  a  function  of  packing  density  and 
extracellular  volume  fraction  (31).  Therefore,  the  total  water 
signal  in  the  intermediate-time  regimen  reflects  a  superposi¬ 
tioning  of  hindered  and  restricted  water,  leading  to  a  non¬ 
monoexponential  signal  decay  at  high  b-  value.  As  we  see  below, 
separating  the  restricted  water  fraction  from  hindered  and  free 
water  fraction  provides  a  more  sensitive  and  specific  biomark¬ 
er  for  tumor  cellularity  compared  with  traditional  ADC. 

Diffusion  models 

Diffusion  models  form  the  fundamental  basis  through  which 
quantitative  information  about  the  underlying  tissue  micro¬ 
structure  can  be  gleaned  from  DWI  signals.  The  typical  milli- 
metric  scale  at  which  DWI  measurements  are  made  (i.e.,  voxel 
size)  is  large  compared  with  the  micrometric  scale  of  the 
underlying  physical  diffusion  process,  and,  therefore,  the  goal 
of  the  diffusion  model  is  to  bridge  the  gap  such  that  inference 
can  be  drawn  on  a  scale  much  smaller  than  the  voxel  dimen¬ 
sions.  One  of  the  advantages  of  DWI  over  other  MRI  techniques 
is  that  the  physical  scale  probed  by  the  measurements  can  be 
adjusted  by  the  experimental  variables,  namely  the  diffusion 
time  (7d)  and  diffusion- weighting  factor  (fo-value).  As  we  see, 
this  forms  the  fundamental  basis  through  which  more 
advanced  multiscale,  or  non-Gaussian,  diffusion  models  offer 
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insight  into  compartmental  diffusion  in  cancer  and  improved 
imaging  biomarkers  for  tumors.  We  begin  this  section  with 
brief  introduction  of  the  classic  ADC  and  tensor  model  for 
Gaussian  diffusion  and  follow  it  with  a  more  detailed  intro¬ 
duction  to  RSI,  a  technique  developed  in  our  laboratory  for 
probing  tissue  microstructure  in  greater  detail. 

ADC.  As  a  departure  from  earlier  studies  in  which  efforts 
were  made  to  measure  the  true  diffusion  process  in  biologic 
systems  (37),  it  was  suggested  in  the  mid  1980s  (24)  to  model 
the  complex  diffusion  in  tissue  using  the  free  (Gaussian) 
diffusion  equation,  but  replace  the  intrinsic  diffusion  coefficient 
D  with  a  global  statistical  parameter  called  ADC.  The  ADC 
concept  has  since  been  used  extensively  in  the  literature  as  well 
as  clinical  oncology  as  a  surrogate  marker  of  tumor  cellularity 
(1-3).  In  practice,  the  ADC  is  estimated  by  combining  the 
experimental  variables  (i.e„  the  magnitude,  duration,  and  tem¬ 
poral  spacing  of  the  diffusion  gradients)  into  a  single  parameter 
called  the  diffusion-weighting  factor,  or  6-value  (24),  and  com¬ 
paring  the  signal  attenuation  at  one  or  more  nonzero  6-values 
with  the  baseline  signal  measured  without  diffusion  weighting 
(i.e.,  with  a  6-value  of  zero).  For  Gaussian  diffusion,  the  signal 
attenuation  decays  exponentially  with  the  product  of  the  6- 
value  and  ADC,  and,  therefore,  the  ADC  reflects  the  slope  of  the 
best  fit  line  to  the  log  signal  as  a  function  of  6-value.  The 
diffusion  tensor  model  (DTI)  extends  the  ADC  concept  to  three- 
dimensional  space,  allowing  for  different  ADCs  along  and 
perpendicular  to  the  principal  axis  of  diffusion  (38). 

The  main  limitation  of  the  ADC  and  tensor  model  as  imaging 
biomarkers  for  tumors  is  their  strict  dependence  on  a  single 
Gaussian  function  for  the  displacement  distribution  of  water 
molecules  within  cancer  tissue  and,  therefore,  a  monoexpo¬ 
nential  dependence  on  the  6-value.  Although  a  single  Gaussian 
assumption  may  be  appropriate  for  data  collected  over  a  fairly 
narrow  range  of  6-values  (up  to  about  1,000  s/mm2),  the 
inferences  that  can  be  drawn  are  limited  by  lack  of  specificity. 
For  example,  numerous  investigators  have  shown  a  negative 
correlation  between  the  tumor  ADC  values  and  cellularity 
(1-3),  which  is  commonly  attributed  to  increased  restricted 
diffusion  imposed  by  tumor  cells.  However,  despite  increased 
restricted  diffusion,  tumor  ADC  values  rarely  fall  below  that  of 
normal  appearing  white  matter  (NAWM).  This  is  true  even  in 
highly  cellular  tumors  that  originate  in  white  matter,  such  as 
glioblastoma  multiforme  (GBM)  and  primary  CNS  lymphoma 
(39).  One  explanation  for  higher  than  expected  ADC  values  in 
these  tumors  is  the  presence  of  vasogenic  edema  and  focal 
necrosis  within  the  tumor  itself,  which  increases  the  ADC 
through  reduced  hindrance  imposed  on  the  extracellular  water 
(40,  41).  Thus,  increased  ADC  due  to  edema  and  necrosis  will 
offset  reduced  ADC  imposed  by  tumor  cells,  resulting  in  lesions 
that  are  difficult  to  discern  from  NAWM  on  the  basis  of  the 
ADC  alone.  This  effect  becomes  more  pronounced  as  the  6- 
value  is  reduced  because  of  increasing  sensitivity  to  the  fast, 
hindered  water  fraction. 

RSI.  Over  the  past  decade  or  so,  advances  in  DWI  acqui¬ 
sition  and  gradient  hardware  have  made  it  possible  to  probe  a 
whole  new  regimen  of  water  diffusion  in  cancer  tissue  beyond 
what  was  previously  possible  on  clinical  MRI  scanners.  Spe¬ 
cifically,  the  application  of  strong  magnetic  field  gradient 


pulses  has  allowed  a  much  greater  range  of  6-values,  diffusion 
directions,  and  diffusion  times  to  be  acquired  during  a  clinical 
acquisition.  As  a  result,  diffusion  signals  can  be  made  specific 
to  pools  of  water  with  very  low  effective  diffusion  coefficients, 
well  below  the  tortuosity  limit  for  ECS  water,  likely  originating 
from  restricted  water  trapped  within  the  cells  themselves  (42). 
Numerous  promising  methods  are  emerging  to  capture  and 
model  complex  non-Gaussian  diffusion  in  tissue,  including 
biexponential  (20, 21),  stretched  exponential  (22),  and  kurtosis 
models  (23),  as  well  as  methods  to  study  perfusion-related 
effects  in  DWI  data  at  low  6-values,  such  as  IVIM  (24).  The 
application  and  description  of  many  of  these  techniques  for 
neurooncology  applications  can  be  found  in  an  excellent 
review  provided  by  Maier  and  colleagues  (39).  Moreover, 
techniques  are  emerging  to  probe  diffusion  spectra  at  ultra- 
short  diffusion  times  on  clinical  systems  using  oscillating 
gradients  (43),  which  hold  promise  for  probing  intracelluar 
structures  and  alterations  in  cancer  (44).  For  this  review,  we 
focus  on  one  particular  technique  developed  in  our  laboratory 
called  RSI.  Over  the  last  few  years,  RSI  has  gained  increased 
recognition  as  an  important  tool  in  oncology  that  overcomes 
many  of  the  limitation  of  traditional  DWI  and  ADC. 

RSI  is  a  general  framework  for  modeling  diffusion  signals 
collected  across  a  broad  range  of  experimental  parameters  and 
relating  these  signals  to  underlying  tissue  parameters  (e.g.,  size 
and  shape/ orientation  of  hindered  and  restricted  water  com¬ 
partments)  using  a  linear  mixture  model  (45).  The  ultimate 
goal  of  RSI  is  to  enable  quantitative  estimates  of  tissue  micro¬ 
structure  based  on  noninvasive  imaging.  To  achieve  this,  the 
diffusion  signal  is  modeled  as  reflecting  a  mixture  of  compo¬ 
nents,  where  each  component  describes  the  signal  dependence 
on  specific  tissue  properties  (e.g.,  cell  size,  density,  orientation, 
etc.)  as  a  function  of  the  experimental  (protocol)  settings  (e.g., 
6-value,  diffusion  time,  echo  time,  etc).  The  total  signal 
becomes  the  weighted  sum  of  these  components,  and  the  goal 
is  to  determine  the  individual  weights.  This  is  achieved  through 
the  application  of  generalized  linear  estimation  techniques 
(45,  46).  As  such,  the  RSI  framework  is  designed  to  strike  a 
balance  between  model  complexity  and  interpretability  by 
minimizing  a  priori  assumptions  on  microstructure  while  pre¬ 
serving  biophysical  interpretability  of  the  resulting  estimates. 

Our  current  clinical  implementation  of  RSI  acquires  data 
with  6-values  of  500,  1,500,  and  4,000  s/mm2  and  multiple 
diffusion  directions  at  each  6-value  at  a  fixed  intermediate 
diffusion  time  (~90  ms),  where  intracellular  spins  would  be 
expected  to  be  in  the  fill-up  regimen.  The  RSI  design  matrix 
includes  a  distribution  (or  "spectrum")  of  effective  diffusion 
pools  spanning  hindered  and  restricted  length  scales  with  both 
isotropic  and  anisotropic  geometries.  In  Fig.  3,  we  illustrate 
both  a  schematic  for  the  RSI  spectrum  model  and  the  resultant 
fit  of  the  model  to  data  collected  in  a  51  year-old  patient  with 
right  frontal  GBM  before  surgical  intervention  (47).  Note  the 
separation  of  diffusion  components  in  different  tissue  types, 
with  the  primary  lesion  exhibiting  a  large  volume  fraction  of 
spherically  restricted  water,  likely  stemming  from  water 
trapped  within  cancer  cells.  Also  note  the  large  fraction  of 
free  and  hindered  water  in  areas  of  necrosis  and  edema.  The 
volume  fraction  of  spherically  restricted  water  (either  with  or 
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Figure  3.  RSI  analysis  of  a  51  year 
old  male  with  right  frontal  GBM.  A, 
illustration  of  the  RSI  "spectrum" 
model  used  to  fit  the  multi-b-value, 
multidirection  DWI  data.  Scales  0-2 
and  3-6  correspond  to  restricted 
and  hindered  diffusion, 
respectively.  Scales  0,  6,  and  7  are 
isotropic,  whereas  scales  1-5  are 
anisotropic  (i.e.,  oriented). 

B,  RSI-derived  OV  weighted) 
volume  fraction  maps  for  each 
scale  in  A.  C,  Ti -weighted 
postcontrast  (D)  T2-weighted 
FLAIR  (E)  RSI-derived  "cellularity 
map"  (RSI-CM)  corresponding  to  a 
weighted  ("beamformed")  linear 
combination  of  scales  0-7  showing 
maximal  sensitivity  and  specificity 
to  spherically  restricted  diffusion 
(scale  0).  F,  bar  plot  of  volume 
fractions  for  two  representative 
voxels  in  tumor  and  necrotic  tissue, 
respectively. 
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without  additional  filtering  of  white  matter  signal  using  a 
technique  called  "beamforming";  ref.  48)  has  been  coined  the 
RSI  cellularity  index  or  cellularity  map  (47, 49-51).  In  addition, 
the  signal  from  the  restricted  anisotropic  component  can  be 
used  to  map  the  density  and  orientation  of  surrounding  white 
matter  tracts  (50). 

In  summary,  advanced  DWI  methods  such  as  RSI  that 
acquire  data  over  an  extended  i-value  range  provide  the  ability 
to  quantify  complex  non-Gaussian  diffusion  in  tissue  (21- 
23,  45,  52).  Although  these  emerging  techniques  offer  a  new 
class  of  cancer  imaging  biomarkers,  there  remain  many  unan¬ 
swered  questions  and  tremendous  opportunity  for  further 
advancing  the  field.  Of  particular  importance  is  to  understand 
exactly  how  parameters  of  the  diffusion  model  relate  to  specific 
properties  of  cancer  tissue,  such  as  tumor  cell  size,  density,  and 
nuclear  volume  fraction.  How  specific  are  these  biomarkers  to 
tumor  cells  versus  healthy  or  inflammatory  cells?  Can  these 
new  imaging  biomarkers  characterize  tumor  cytotypes  and 


address  a  top  challenge  posed  by  the  NCI  to  develop  in  vivo 
imaging  methods  that  can  probe  the  identity,  quantity,  and 
location  of  different  cells  that  make  up  a  tumor  and  its 
microenvironment?  (NIH  RFA-CA-13-020).  Finally,  to  what 
extent  can  these  new  imaging  biomarkers  help  solve  current 
clinical  dilemmas  in  oncology  such  as  distinguishing  tumor 
recurrence  from  radiation  injury?  Answering  these  questions 
will  undoubtedly  require  both  advances  in  diffusion  modeling, 
simulation,  and  MRI  hardware  and  software  [such  as  the  latest 
ultra-high  performance  ("connectome")  MRI  scanners;  ref.  53] 
together  with  improved  quantitative  histology,  multimodal 
image  registration,  and  validation  procedures. 

Clinical  Applications  of  RSI 

Because  of  its  ability  to  isolate  areas  of  truly  restricted 
diffusion  by  separating  and  removing  the  hindered  diffusion 
signal,  RSI  offers  a  more  direct  measure  of  tumor  cells  than 
other  diffusion-weighted  methods.  Within  the  past  year,  we 
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Figure  4.  From  left  to  right,  Ti- 
weighted  postcontrast,  T2- 
weighted  FLAIR,  ADC,  and  RSI 
cellularity  map  for  a  53  year  old 
male  with  treatment-naive  right 
temporal  GBM  (top)  and  a  73  year- 
old  female  with  metatstatic  non- 
small-cell  lung  cancer.  Bottom, 
ROC  curves  demonstrating 
increased  sensitivity,  specificity, 
and  overall  accuracy  for  delineating 
high-grade  primary  and  metastatic 
brain  tumors  with  RSI  compared 
with  ADC.  Note  the  high  tumor 
conspicuity  on  RSI  and  the  more 
protruding  finger-like  margins  in 
GBM  compared  with  metastatic 
disease,  consistent  with  infiltrating 
tumor  into  peritumoral  edema. 


have  shown  that  by  identifying  areas  of  restricted  diffusion,  RSI 
can  be  applied  clinically  to  (i)  improve  tumor  conspicuity  in 
patients  with  high-grade  primary  and  metastatic  brain  tumors 
(47),  (ii)  facilitate  the  interpretation  of  tumor  response  on 
imaging  following  antiangiogenic  treatments  that  significantly 
reduce  edema  (49),  and  (iii)  improve  the  ability  to  visualize 
white  matter  pathways  coursing  through  regions  of  peritu¬ 
moral  edema  relative  to  standard  DTI  (50).  In  the  following 
sections,  we  describe  each  of  these  clinical  applications  in 
further  detail.  In  addition,  we  describe  a  clinical  scenario  in 
which  RSI  lacks  specificity  and  a  multispectral  imaging 
approach  is  warranted. 


Improved  conspicuity  of  high-grade  tumors  with  RSI 

As  described  in  the  Introduction,  ADC  is  frequently  used  as  a 
marker  of  tumor  cellularity  in  patients  with  high-grade  tumors 
(1-3).  Areas  of  tumor  are  associated  with  decreased  ADC 
relative  to  surrounding  tissue  (40).  However,  concomitant 
edema  and  tumor-related  necrosis  increase  ADC  values,  there¬ 
by  directly  opposing  the  reduction  in  ADC  associated  with 
tumor  (54,  55).  This  offset  presents  a  diagnostic  challenge  by 
diminishing  the  conspicuity  of  tumor  on  ADC  maps.  Because 
RSI  isolates  areas  of  spherically  restricted  diffusion,  we  tested 
whether  it  could  provide  increased  conspicuity  and  delineation 
of  tumor  margins  relative  to  standard  and  high  ft- value  ADC 
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(47).  To  accomplish  this,  RSI  was  performed  in  ten  presurgical 
patients:  four  with  GBM,  three  with  primary  CNS  lymphoma, 
and  three  with  metastatic  brain  tumors.  Tumor  conspicuity, 
edema  conspicuity,  and  relative  sensitivity  to  edema  were 
quantified  for  RSI  cellularity  maps  (RSI-CM),  high  lvalue  DWI 
(b  =  4,000),  and  ADC,  and  these  values  were  compared  in 
manually  drawn  volumes  of  interest.  Receiver  operating  char¬ 
acteristic  (ROC)  curves  were  used  to  evaluate  the  sensitivity 
and  specificity  of  each  method  for  delineating  tumor  from 
NAWM.  In  addition  to  visible  differences  in  conspicuity  (Fig.  4), 
ROC  curves  revealed  greater  sensitivity  and  specificity  for 
delineating  tumor  from  NAWM  with  RSI-CM  (AUC  =  0.91) 
compared  with  both  high  6-value  DWI  (AUC  =  0.77)  and  ADC 
(AUC  =  0.66).  In  addition,  the  relative  sensitivity  to  edema  was 
greater  for  high  6-value  DWI  and  ADC  compared  with  RSI, 
reflecting  RSI's  ability  to  suppress  the  fast  diffusion  component 
associated  with  edema.  Furthermore,  greater  heterogeneity  of 
the  diffusion  signal  within  the  tumor  was  observed  on  the  RSI- 
CMs  compared  with  DWI  and  ADC,  as  evidenced  by  a  broader 


histogram  distribution.  This  may  represent  the  intrinsic  het¬ 
erogeneity  of  tumor  cellularity  both  within  and  across  tumor 
types.  These  data  demonstrate  one  promising  application  of 
RSI,  i.e.,  improved  conspicuity  and  delineation  of  high-grade 
tumors  compared  with  traditional  DWI  models  and  under¬ 
score  the  possibility  that  RSI  may  prove  helpful  in  delineating 
tumor  cytotypes  and  infiltrating  disease  in  peritumoral  edema. 

RSI  in  the  context  of  antiangiogenic  treatment 

Antiangiogenic  therapies,  such  as  bevacizumab,  are  increas¬ 
ingly  used  in  the  treatment  of  recurrent  high-grade  gliomas. 
However,  these  agents  decrease  permeability  of  the  blood- 
brain  barrier  and,  therefore,  decrease  contrast  enhancement 
and  edema  in  patients  with  high-grade  gliomas  in  a  manner 
that  may  not  correlate  with  actual  tumor  response — a  phe¬ 
nomenon  known  as  pseudoresponse  (56).  Given  this  imaging 
challenge,  we  evaluated  the  ability  of  RSI  to  improve  conspi¬ 
cuity  within  regions  of  the  tumor  compared  with  ADC  in 
patients  treated  with  bevacizumab  and  to  further  demonstrate 


Figure  5.  A  67  year  old  male  with  left  parietal  GBM  status  postresection  and  chemoradiation.  Top,  theTi  postcontrast  -  Ti  precontrast  (A),  FLAIR  (B),  ADC  (C), 
and  RSI-CMs  (D)  before  the  start  of  bevacizumab;  middle,  Ti  postcontrast  -  Ti  precontrast  (E),  FLAIR  (F),  ADC  (G),  and  RSI-CMs  (H)  after  initiation  of 
bevacizumab.  Arrowheads,  contrast-enhancing  region  (green),  the  surrounding  region  of  FLAIR  hyperintensity  (yellow),  and  the  region  of  restricted  diffusion 
on  RSI-CMs  (red).  Although  there  is  a  decrease  in  contrast  enhancement  and  surrounding  FLAIR  hyperintensity  after  initiation  of  bevacizumab,  the  region  of 
restricted  diffusion  increases  and  becomes  more  confluent,  suggesting  worsening  residual/recurrent  tumor.  Moreover,  this  increase  in  the  region  of  restricted 
diffusion  is  much  more  conspicuous  on  the  RSI-CMs  compared  with  the  ADC.  Bottom  row  depicts  these  changes  on  "change  maps"  (change  in  Ti 
postcontrast  -  precontrast;  I),  change  in  FLAIR  (J),  change  in  ADC  (K),  and  change  in  the  RSI-CMs  (L),  with  red-yellow  indicating  an  increase  in  signal 
intensity  and  blue-cyan  indicating  a  decrease  in  signal  intensity.  Of  note,  on  the  ADC  change  map  (K),  the  area  of  increased  restricted  diffusion  is  essentially 
masked  by  the  decreased  signal  intensity  within  the  region  of  surrounding  FLAIR  hyperintensity. 
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that  RSI  is  minimally  affected  by  bevacizumab-induced  reduc¬ 
tions  in  edema  (49).  RSI-CMs  and  DWI  were  available  for  a 
series  of  patients  with  recurrent  gbomas  at  baseline  and 
following  initiation  of  bevacizumab.  Results  showed  that  all 
patients  exhibited  sharp  decreases  in  contrast  enhancement 
and  edema  following  treatment  (i.e.,  pseudoresponse).  Beva¬ 
cizumab-induced  decreases  in  edema  had  a  greater  effect  on 
ADC  than  on  the  RSI-CMs,  with  the  relative  sensitivity  to 
changes  in  edema  being  more  than  20  times  higher  on  ADC 
than  on  RSI-CMs  (Fig.  5).  These  data  provide  additional 
evidence  that  RSI  is  less  influenced  by  changes  in  edema 
compared  with  ADC,  which  may  confer  an  advantage  of  RSI 
for  interpreting  true  tumor  response  in  the  setting  of  anti- 
angiogenic  treatment.  However,  this  study  did  not  include 
clinical  follow-up.  Thus,  whether  RSI  has  greater  predictive 
validity  than  ADC  is  of  key  importance  and  requires  additional 
investigation. 

RSI  for  improved  delineation  of  white  matter  tracts 

There  is  increasing  enthusiasm  for  the  use  of  DTI  and 
tractography  in  neurosurgical  planning  (57-64),  and  there  is 
some  evidence  that  tractography-guided  neuronavigation  can 
be  used  to  minimize  neurologic  morbidity  (60,  65,  66).  How¬ 
ever,  the  ability  to  resolve  white  matter  structure  in  peritu- 
moral  regions  that  include  edema  has  remained  a  challenge 
(39, 67).  We  applied  the  same  logic  as  in  the  previous  examples 
and  tested  the  ability  of  RSI  to  provide  better  visualization  and 
quantification  of  white  matter  tracts  in  regions  that  include 
edema  (50).  In  this  example,  isolating  the  slow,  restricted 
compartment  yields  a  better  estimate  of  "tubularity"  (i.e., 


models  the  cylindrically  restricted  diffusion  within  axons; 
Fig.  3B,  scale  1;  ref.  45).  This  increases  the  sharpness  of  the 
estimates,  allowing  for  better  delineation  of  fiber  tract  orien¬ 
tation.  In  a  series  of  ten  patients  with  high-grade  gliomas,  we 
were  able  to  demonstrate  that  RSI  yielded  higher  fractional 
anisotropy  (FA)  estimates  in  regions  of  edema  relative  to 
standard  DTI.  Furthermore,  at  follow-up  when  the  edema  had 
resolved  in  most  patients,  FA  estimates  increased  with  DTI,  but 
remained  stable  with  RSI,  indicating  that  FA  estimates  based 
on  DTI  were  artificially  suppressed  by  the  edema.  Tractogra¬ 
phy  performed  within  regions  of  edema  revealed  superior 
ability  of  RSI  to  track  fibers  through  areas  of  significant  edema 
relative  to  standard  DTI  (Fig.  6).  These  data  address  yet  another 
important  and  growing  application  of  diffusion  imaging  within 
the  field  of  neurooncology  (i.e.,  surgical  planning)  and  speak  to 
the  advantage  of  using  advanced  DWI  models  for  revealing  the 
anatomic  structure  of  peritumoral  white  matter.  An  impor¬ 
tant  caveat  to  mention  is  that  it  is  well  known  that  gliomas 
grow  via  an  infiltrative  pattern  and  that  there  is  often 
nonenhancing  infiltrating  tumor  surrounding  the  enhancing 
portion  of  a  high-grade  glioma  (68).  Further  investigation  is 
required  to  determine  whether  RSI  or  other  advanced  DWI 
methods  will  prove  useful  for  detecting  these  areas  of  tumor 
infiltration. 

Importance  of  a  multispectral  imaging  approach 

Despite  the  many  recent  advancements  in  MR  technology 
and  DWI  methods,  no  single  modality  has  emerged  as  the  "holy 
grail"  of  tumor  imaging.  Thus,  most  diagnostic  challenges  in 
neurooncology  are  best  addressed  using  a  multispectral 
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Figure  6.  Streamline  tractography  of  the  superior  longitudinal  fasciculus  for  a  58  year  old  female  with  a  right  temporal  lobe  GBM  projected  onto  baseline  and 
follow-up  FLAIR  images.  Left,  RSI  and  DTI-based  tractography  at  baseline  in  regions  of  edema;  right,  data  obtained  using  the  same  tractography  algorithm 
once  the  edema  had  mostly  resolved.  The  ipsilateral  (red)  and  contralateral  (green)  3D  renditions  of  the  superior  longitudinal  fasciculus  are  superimposed  on 
axial  and  sagittal  FLAIR  slices  collected  at  each  time  point.  The  GBM  is  shown  in  blue  in  the  preoperative  image.  With  RSI,  the  superior  longitudinal 
fasciculus  appears  very  similar  at  baseline  and  at  follow-up.  However,  with  DTI,  the  superior  longitudinal  fasciculus  appears  thinner  and  truncated  at  baseline 
in  regions  of  edema.  Black  arrows,  frontal  and  parietal  regions  of  the  superior  longitudinal  fasciculus  that  terminate  completely  in  regions  of  edema;  red  arrow, 
sparse  streamlines  in  the  temporal  portion  of  the  superior  longitudinal  fasciculus.  These  streamlines  are  "recovered"  using  DTI  once  the  edema  resolves. 
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Figure  7.  Comparison  of  a  55  year 
old  male  with  GBM  treated  with 
chemoradiation  and  bevacizumab 
(top)  and  a  66  year  old  male  with 
GBM  before  any  treatment 
(bottom).  Ti  postcontrast  images 
(A  and  D),  RSI  maps  (B  and  E),  and 
rCBV  maps  (C  and  F)  are  shown. 
Degree  and  homogeneity  of 
restricted  diffusion  is  greater  in  the 
patient  treated  with  bevacizumab 
than  in  the  pretreatment  GBM 
control  (images  scaled  identically 
with  same  window  and  level), 
whereas  rCBV  in  the  region  of 
restricted  diffusion  is  remarkably 
low  in  the  patient  treated  with 
bevacizumab— lower  than  in  the 
GBM  control  and  lower  than  in  the 
NAWM. 


imaging  approach.  One  particular  challenge  occurs  following 
treatment  with  concurrent  bevacizumab  and  radiotherapy, 
with  a  subset  of  patients  developing  regions  of  marked  and 
persistent  restricted  diffusion  that  do  not  seem  to  reflect  an 
aggressive  tumor  (69-71).  Although  the  etiology  of  these  lesions 
remains  uncertain,  pathologic  confirmation  in  several  patients 
has  revealed  atypical  gelatinous  necrosis.  We  have  coined  this 
abnormality  bevacizumab-related  imaging  abnormality  (BRIA) 
and  have  observed  that  the  BRIA  signal  on  RSI  is  quantitatively 
similar  to  that  seen  in  a  tumor  (51).  This  illustrates  a  clinical 
scenario  in  which  RSI  lacks  specificity  and  a  multispectral 
imaging  approach  is  warranted.  Thus,  we  explored  whether 
using  RSI  in  combination  with  perfusion  imaging  could  help  to 
differentiate  BRIA  from  recurrent  tumor.  In  a  series  of  patients, 
we  show  that  these  techniques  are  complementary  in  that  RSI 
is  superior  to  rCBV  for  differentiating  pathology  from  NAWM, 
whereas  rCBV  is  superior  for  differentiating  BRIA  from  tumor 
(Fig.  7).  Thus,  the  combination  of  high  RSI  signal  and  low  rCBV 
provides  a  distinct  imaging  signature  of  the  BRIA  phenomenon. 
Because  of  the  increasing  use  of  antiangiogenic  agents,  imaging 
methods  that  increase  our  understanding  of  both  pseudore¬ 
sponse  and  BRIA  are  of  high  importance.  These  data  highlight 
an  important  example  in  which  RSI  lacks  specificity  in  its 
current  instantiation,  and  information  from  complementary 
imaging  modalities  is  essential. 

RSI  for  improved  detection  of  other  solid  organ  tumors 

Although  this  review  is  primarily  focused  on  neurooncology 
applications,  quantitative  imaging  is  equally  relevant  to  other 


solid  organ  tumors.  For  example,  multiparametric  MRI  has 
been  explored  in  a  variety  of  applications,  including  discrim¬ 
ination  between  indolent  and  aggressive  disease  in  prostate 
cancer.  Routine  r2-weighted  imaging  of  the  prostate  is  the 
most  sensitive  way  to  evaluate  anatomic  detail,  but  is  limited 
for  disease  detection,  with  sensitivity  of  around  70%  and 
specificity  of  55%  (72).  Functional  MR  techniques  enhance 
detection,  grading,  and  staging  of  prostate  cancer  through  the 
use  of  dynamic  contrast  enhancement  (DCE),  DWI,  and  MR 
spectroscopic  imaging.  DCE  requires  intravenous  administra¬ 
tion  of  a  T ,-shorlenmg  agent  with  the  chief  limitation  being 
that  the  maximum  contrast  enhancement  between  malignant 
and  nonaggressive  disease  is  small  (73).  MR  spectroscopic 
imaging  increases  specificity  (74,  75),  but  is  technically  chal¬ 
lenging  and  can  add  significantly  to  scan  time. 

Multiple  studies  have  shown  that  DWI  improves  sensitivity 
and  specificity  in  the  diagnosis  of  prostate  cancer  by  increasing 
tumor  conspicuity  on  DWI  or  quantitative  ADC  maps.  How¬ 
ever,  hemorrhage,  inflammatory  processes,  and  benign 
nodules  in  the  transitional  zone  can  all  exhibit  lower  ADC 
values,  leading  to  false  positives  (76).  DWI  can  also  suffer  from 
severe  spatial  distortion,  limiting  its  coregistration  to  anatom¬ 
ic  images,  which  is  necessary  for  tumor  localization. 

Increasing  Gleason  score  correlates  with  loss  of  normal 
gland  formation,  loss  of  peripheral  gland  tubular  structure, 
and  increased  cellularity  (77).  We  hypothesize,  therefore,  that 
RSI  cellularity  will  correlate  with  higher  tumor  grade,  as 
measured  by  Gleason  score,  and  will  provide  significantly 
greater  accuracy  in  discriminating  aggressive  tumors  from 
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Figure  8.  Gleason  3+4.  A,  histology 
section  stained  with  hematoxylin 
and  eosin.  Blue  dotted  line,  the 
boundary  of  the  tumor.  B,  RSI 
cellularity  map,  color-coded  and 
overlaid  onT2.  C,  ADC  image.  D,T2 
image.  E,  3D  volume  rendering  of 
the  RSI  (in  yellow),  the  whole  extent 
of  the  prostate  as  traced  on  T2 
images  (translucent  blue),  and 
green  lines  indicating  the  boundary 
of  the  tumor  on  each  of  the  whole- 
mount  histologic  sections  that 
were  compared  with  the  RSI.  The 
white  arrow  indicates  the  line 
corresponding  to  the  histology 
section  shown  in  A.  F,  raw 
perfusion  data. 


benign  and  indolent  lesions  when  compared  with  current 
functional  or  anatomic  imaging  techniques. 

Our  preliminary  data  are  encouraging.  Figure  8  shows  a 
patient  with  Gleason  3+4  =  7  disease.  RSI  cellularity  map  is 
shown  color-coded  and  fused  with  the  T2  after  spatial  distor¬ 
tion  correction  and  registration.  Note  the  correspondence  of 
RSI  cellularity  with  the  histopathology  slice.  T2  and  perfusion 
images  are  not  as  conspicuous,  with  much  less  contrast-to- 
noise  ratio.  This  represents  an  example  case  in  which  conven¬ 
tional  imaging  is  not  as  diagnostic. 

Importance  of  spatial  distortion  correction  for  accurate 
image-guided  intervention 

One  of  the  main  limitations  of  diffusion  imaging  in  general  is 
nonlinear  spatial  distortion  of  the  images  due  to  a  number  of 
factors,  including  gradient  nonlinearities,  eddy  currents,  and  B0 
held  inhomogeneities.  Although  distortions  due  to  nonlinea¬ 
rities  of  the  gradient  Helds  are  commonly  corrected  for  by 
software  on  the  scanner  console,  B0  Held  inhomogeneities  are 
not,  despite  being  the  dominant  source  of  spatial  inaccuracy  in 
DWI.  The  magnitude  of  B0  distortions  varies  depending  on  a 
number  of  factors,  including  field  strength,  positioning  of  the 
subject  within  the  scanner,  and  subject-specific  anatomy.  The 
typical  pattern  of  B0  distortion  on  3T  systems  is  illustrated 


in  Fig.  9,  along  with  the  total  whole-brain  histogram  of  root 
mean  square  (RMS)  displacements.  As  illustrated  in  the  figure, 
on  average  (based  scans  from  on  40  subjects)  the  typical 
(mode)  distortion  magnitude  is  approximately  2  mm,  with  a 
substantial  proportion  of  voxels  displaced  more  than  6  mm. 
Such  distortions  are  of  particular  concern  when  images  are 
used  for  image-guided  intervention,  including  surgery,  biopsy, 
or  radiation  dose  planning.  Although  methods  for  correcting  B0 
distortions  have  existed  for  some  time  (78-80),  they  are  typ¬ 
ically  used  only  in  research  studies  and  not  in  clinical  practice, 
primarily  due  to  the  additional  scan  time  required  to  acquire 
the  B0  field  maps  required  for  standard  correction  methods. 
The  method  introduced  by  Holland  and  colleagues  (81)  over¬ 
comes  this  limitation  by  requiring  only  a  single  additional  TR 
(2-3  seconds),  using  the  reverse  phase-encode  polarity  meth¬ 
od  (79,  82,  83).  Moreover,  this  additional  volume  acquisition 
can  be  integrated  directly  in  the  native  DWI  protocol  without 
changing  the  pulse  sequence.  Combined  with  efficient  post¬ 
processing  methods  to  estimate  the  distortion  field  (81),  the 
reverse  phase-encode  polarity  technique  provides  a  much 
needed  clinical  solution  for  accurate  spatial  distortion  correc¬ 
tion  of  DWI  data.  The  RSI  technique  incorporates  this  proce¬ 
dure  as  part  of  standard  preprocessing  of  data,  and  all  the  data 
presented  in  this  review  has  been  corrected  in  this  manner. 
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Figure  9.  Three-plane  map  and 
histogram  of  the  RMS  displacement 
of  voxels  (in  mm)  due  to  Bo 
distortions  for  axially  acquired  EPI 
images  (A/P  phase  encoding)  from 
EPI  scans  of  a  40  normal  healthy 
subjects. 
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Discussion  and  Conclusion 

The  field  of  MR  imaging  is  rapidly  evolving,  leading  to  new 
and  exciting  possibilities  within  neurooncology,  urologic 
oncology,  and  beyond.  In  particular,  the  development  of 
advanced  DWI  methods  has  allowed  for  improved  visualiza¬ 
tion  and  detection  of  tumor  cells  and,  thus,  has  great  potential 
for  better  understanding  of  tumor  biology.  Numerous  clinical 
applications  of  this  powerful  technique  have  already  been 
demonstrated,  including  tumor  characterization  and  grading, 
prognostication,  early  prediction  of  response  to  therapy  and 
survival,  distinguishing  tumor  from  treatment-related  con¬ 
founds,  detecting  microinfiltration,  and  guiding  neurosurgical 
and  radiation  planning  (84).  However,  future  advances  in  the 
field  will  require  a  fundamental  understanding  of  the  under¬ 
lying  DWI  signal  coupled  with  validation  of  diffusion  contrast 
in  unique  tumor-related  pathologies.  In  addition  to  improved 
detection  of  tumors,  advanced  diffusion  methods  such  as  RSI 
may  also  provide  quantitative  characterization  of  cellular 
properties  such  as  cell  size,  permeability,  and  nuclear  volume 
fraction,  based  on  signal  variation  as  a  function  of  diffusion 


time  and  echo  time  (42).  Validation  of  these  measures  will 
likely  be  borne  out  of  translational  efforts  that  include  both 
preclinical  and  clinical  studies  in  which  histologic  specimens 
are  carefully  coregistered  to  in  vivo  imaging.  The  need  for 
targeted  biopsies  based  on  advanced  DWI,  precise  coregistra¬ 
tion  of  DWI  with  other  imaging  modalities,  and  careful  correc¬ 
tions  for  geometric  distortions  will  all  be  pivotal  to  providing 
the  spatial  precision  needed  to  achieve  such  validation.  These 
requirements  underscore  the  need  for  a  multidisciplinary 
approach  to  this  validation  including  experts  in  the  fields  of 
oncology,  surgery,  pathology,  and  radiology.  Through  this  col¬ 
laboration,  not  only  will  current  applications  of  DWI  be  further 
improved,  but  new  possibilities  will  also  be  created  that  will 
ultimately  lead  to  better  care  for  patients  suffering  from  cancer. 
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Abstract 


Background:  Standard  magnetic  resonance  imaging  (MRI)  of  the  prostate  lacks  sensitivity  in  the 
diagnosis  and  staging  of  prostate  cancer  (PCa).  To  improve  the  operating  characteristics  of 
prostate  MRI  in  the  detection  and  characterization  of  PCa,  we  developed  a  novel,  enhanced  MRI 
diffusion  technique  using  restriction  spectrum  imaging  (RSI-MRI). 

Methods:  We  compared  the  efficacy  of  our  novel  RSI-MRI  technique  to  standard  MRI  for 
detecting  extra-prostatic  extension  in  7  PCa  patients  with  histologically  proven  pT3a  disease 
from  a  total  of  28  PCa  patients  imaged.  All  28  patients  underwent  MRI  and  RSI-MRI  prior  to 
radical  prostatectomy. 

Results:  Pre-operative  standard  MRI  correctly  identified  extra-prostatic  extension  in  2  of  7  (29%) 
patients,  while  RSI-MRI  identified  extra-prostatic  extension  in  7  of  7  (100%)  patients. 
Conclusions:  We  conclude  that  our  novel  RSI-MRI  technology  shows  promise  for  substantially 
improving  PCa  imaging.  Further  translational  studies  of  prostate  RSI-MRI  in  the  diagnosis  and 
staging  of  prostate  cancer  are  indicated. 


Introduction 


Multiparametric  magnetic  resonance  imaging  (MRI)  is  an  emerging  diagnostic  tool  for  the 
screening,  staging,  and  treatment  of  prostate  cancer  (PCa).  (1-17)  However,  prostate  MRI 
demonstrates  variable  sensitivity  (49  -  88%)  and  specificity  (84-  89%),  which  currently  limits  its 
clinical  utility (2, 3, 14). 

Standard  diffusion  weighted  imaging  (DWI)  improves  the  operating  characteristics  of  prostate 
MRI  (2,3,14,16,18,19).  DWI  detects  the  diffusivity  of  microscopic  water  and  subsequently 
reflects  the  cellularity  and  integrity  of  cells  imaged.  When  combined  with  T2  weighted  imaging 
and  dynamic  contrast  enhanced  imaging,  DWI  improves  sensitivity  and  specificity  of  PCa 
diagnosis  by  increasing  tumor  conspicuity(2,14,16,18,19). 

DWI  also  suffers  from  distortion  due  to  magnetic  field  inhomogeneity  (20)  and  relatively  high 
false  positive  rates  owing  to  hemorrhage  (i.e.  from  prior  biopsy),  inflammatory  processes,  and 
benign  nodules  in  the  transitional  zone  (21).  In  order  to  overcome  these  challenges  and  improve 
magnetic  resonance  detection  of  PCa,  we  have  developed  a  novel,  sophisticated  diffusion  method 
tenned  prostate  Restriction  Spectrum  Imaging  (RSI),  or  RSI-MRI  (22,23).  By  collecting  a 
broader,  more  extended  spectrum  of  diffusion  images(22),  combined  with  sophisticated 
modeling  of  differential  water  compartments  in  tissue  and  correction  of  spatial  distortion(24), 
RSI-MRI  theoretically  focuses  on  the  signal  emanating  from  the  intracellular  water  compartment 
of  tumor  cells  and  thereby  minimizes  false  positive  signals(25).  A  similar  RSI  technique, 


recently  applied  to  brain  imaging,  enhanced  the  signal  of  glioblastoma  multiforme  tumors  by  10- 
fold  (23). 

In  this  pilot  clinical  study,  we  detennined  the  clinical  utility  of  our  novel  prostate  RSI-MRI 
technique  to  stage  PCa  by  assessing  the  accuracy  of  RSI-MRI  to  detect  extra-prostatic  extension 
(EPE)  of  tumor. 

Materials  and  Methods 

RSI-MRI 

In  this  IRB  approved  prospective  study,  28  enrolled  patients  underwent  standard  T2,  perfusion 
(with  Gadolinium)  and  diffusion  protocols  at  3T  (Siemens,  Erlangen  Germany)  with  an 
endorectal  coil  (Table  1).  The  28  patients  had  biopsy  proven  prostate  cancer  and  an  MRI  was 
being  perfonned  for  pre-surgical  planning.  The  restriction  spectrum  protocol  parameters 
included  b-values  0,  800,  1500,  4000  s/mm2  in  30  unique  diffusion  directions  for  each  nonzero  b- 
value.  RSI  cellularity  maps  isolating  the  restricted  isotropic  water  fraction23  were  reconstructed 
based  on  all  b-values  and  then  standardized  across  the  sample  with  z-score  maps.  The  z-score 
maps  were  calculated  by  (1)  measuring  the  mean  and  standard  deviation  of  normal  prostate 
signal  from  the  raw  RSI  cellularity  maps  across  the  patient  pool,  (2)  subtracting  this  measured 
mean  value  from  each  subjects  RSI  cellularity  map  (3)  and  dividing  by  the  measured  normal 
prostate  standard  deviation.  Correction  of  image  distortions  in  RSI  maps  due  to  B0  field 
inhomogeneities  was  performed  using  the  alternating  phase-encode  technique(24).  Per  standard 
clinical  protocol,  radiology  reports  were  generated  by  an  experienced  radiologist  for  each  of  the 
patients  blinded  to  the  RSI  results  (Table  2). 


Correlation  of  RSI-MRI  with  histopathology 

Whole  mount  histopathology  was  performed  on  the  28  enrolled  patients  who  underwent  RSI- 
MRI  and  subsequent  radical  prostatectomy.  For  the  8  radical  prostatectomy  patients  with 
histologically  identified  EPE  (pT3a  disease),  presurgical  pelvic  RSI-MRI  imaging  was  reviewed. 
Out  of  the  8  patients  with  histologically  identified  EPE,  one  patient  was  excluded  due  to  an  error 
in  data  collection  for  the  MRI-RSI,  where  distortion  correction  was  unable  to  be  performed. 
Thus,  a  total  of  7  patients  were  used  for  evaluation  in  this  pilot  study. 

The  area  of  tumor  was  identified  by  an  experienced  uropathologist  (Table  3).  Clinical  images 
were  reviewed  in  a  non-blinded  fashion  to  verify  perfonnance  of  RSI-MRI  (Figure  1),  where 
pathology  was  available  for  comparison.  Standard  MRI  (including  T2,  perfusion,  and  traditional 
DWI),  RSI-MRI,  and  histopathology  were  compared  at  a  multidisciplinary  tumor  board 
consisting  of  radiologists,  urologists,  and  uropathologists.  RSI-MRI  images  were  reviewed  to 
detennine  if  the  RSI  signal  extended  beyond  the  prostatic  capsule,  defined  as  the  T2  hypointense 
thin  line  separating  parenchymal  tissue  from  extraprostatic  fat  and  neurovascular  structures.  RSI 
maps  were  spatially  corrected  and  overlayed  on  standard  T2  maps  to  determine  if  the  RSI  signal 
was  intracapsular  or  bulging,  blurring,  or  grossly  extending  beyond  the  capsule.  This  evaluation 
was  then  compared  to  the  preoperative  MRI  Radiology  reports  to  determine  if  EPE  was 
suspected  without  RSI-MRI. 

Results 


Standard  MRI 


Using  the  standard  preoperative  MR  protocol  (T1  perfusion,  T2,  and  basic  diffusion-weighting), 
patients  1  and  5  demonstrated  radiologic  findings  consistent  with  capsular  involvement  and  EPE. 
Patient  4’s  imaging  findings  were  suspicious  for  EPE  but  not  definitive.  Patients  3,  6,  and  7 
demonstrated  standard  MR  radiologic  findings  consistent  with  bulging  of  the  capsule,  but  no 
indication  of  EPE.  Patient  2’s  standard  imaging  did  not  indicate  concern  for  pathology  abutting 
the  capsule. 

RSI-MRI 

On  analysis  of  RSI-MRI,  7  of  7  (100%)  demonstrated  EPE  (Figure  1,  column  4).  The  areas  of 
EPE  correlated  with  the  areas  encircled  on  the  whole  mount  histopathology.  Four  of  7  (60%) 
patients  with  RSI-MRI  detected  EPE  also  had  positive  surgical  margins.  Figure  2  demonstrates  3 
patients,  out  of  the  total  28  patients  in  the  study,  who  did  not  demonstrate  EPE. 

RSI-MRI  and  Gleason  grading 

As  a  proof  of  principle  in  this  small  sample  of  patients,  regions  of  interest  were  drawn  over  the 
areas  of  the  tumor  on  the  ADC  maps  and  the  RSI-MRI  z-score  maps  to  give  representative  mean 
values  (Table  4).  By  convention,  lower  ADCs  correlate  to  higher-grade  tumors. 

Discussion 

In  this  pilot  clinical  study,  we  demonstrate  the  potential  for  our  novel  prostate  RSI-MRI 
technique  to  substantially  improve  upon  current  imaging-based  modalities  for  the  diagnosis  and 
staging  of  PCa.  In  our  series,  RSI-MRI  successfully  identified  pT3a  disease  in  7  (100%)  of  7 


patients,  while  standard  MRI  accurately  identified  only  2  (29%)  of  7.  These  data  justify  further 
translational  studies  of  RSI-MRI  for  prostate  cancer  detection  and  staging. 

RSI-MRI’s  sensitivity  to  intracellular  diffusivity  differentiates  it  from  standard  DWI(22). 
Extracellular  signal  emanating  from  surrounding  inflammatory  processes  and  hemorrhage 
contribute  less  to  the  signal  detected  by  RSI-MRI  as  compared  to  standard  DWI.  RSI-MRI 
minimizes  these  effects  by  focusing  on  the  signal  emanating  from  intracellular  tumor  cells 
(restricted  diffusion  within  small  spherical  compartments),  with  less  focus  on  the  extracellular 
signal  (hindered  and  free  water  diffusion).  RSI-MRI  achieves  this  by  collecting  diffusion  images 
over  an  extended  b-value  range  (b  =  0  to  4000  s/mm2)  and  decomposing  the  relative  contribution 
from  separable  water  compartments  within  voxels  using  a  linear  mixture  model  framework 
(22,23).  Standard  quantitative  diffusion  techniques  acquire  data  with  b-values  less  than  1,000 
s/mm2.  In  addition,  unlike  conventional  diffusion  imaging  and  ADC  maps,  RSI-MRI  maps  are 
distortion-corrected  and  directly  overlayed  on  T2  anatomic  images  with  voxel  specific  accuracy, 
pennitting  more  accurate  assessment  of  signal  beyond  the  capsular  border.  This  increases  the 
conspicuity  of  the  extraprostatic  extension. 

Current  standard  MRI  prostate  protocols  lack  accurate  localization  of  prostate  cancer.  RSI-MRI 
notably  improves  tumor  localization  with  MRI,  potentially  rendering  MR  imaging  more  relevant 
in  a  variety  of  clinical  scenarios.  Pre-surgical  MRI  results  may  provide  additional  infonnation  for 
clinicians  and  patients  by  informing  surgical  planning  (i.e.  nerve-sparing).  Other  potential 
applications  for  RSI-MRI  meriting  further  study  include  lesion  localization  for  targeting  of 


image-guided  biopsies,  serial  imaging  in  an  active  surveillance  population,  and  evaluation  of 
post-treatment  recurrence  following  primary  radiotherapy. 

A  challenge  in  further  developing  clinical  applications  of  RSI-MRI  lies  in  accurately  defining  its 
specificity,  sensitivity,  and  positive  and  negative  predictive  values.  Limitations  of  this  current 
study  included  the  small  number  of  patients,  and  the  retrospective  non-blinded  nature  of  the 
study.  This  was  necessary  for  this  pilot  study  for  innovative  methodological  development,  but 
suffers  due  to  potential  reader  bias.  An  ongoing,  blinded  prospective  study  will  allow  for  more 
detailed  evaluations  of  the  operating  characteristics  of  RSI-MRI  relative  to  standard  MRI.  RSI- 
MRI  is  a  standardized  quantitative  technique  and  can  potentially  be  employed  across  platforms 
and  institutions.  The  color  bar  in  Figure  1  illustrates  the  z-score  scale  for  the  RSI  maps.  As  the 
number  of  enrollees  increase,  the  z-score  threshold  range  for  concern  for  malignancy  can  be 
detennined  by  pooling  the  variance  across  patients.  One  of  the  challenges  of  standard  diffusion 
imaging,  is  that  the  ADC  values  are  not  standardized  across  MRI  scanners.  The  z-score  is  a 
standardized  statistical  method  and  inherently  normalizes  across  the  patient  pool.  Thus,  the  z- 
score  is  a  value  that  can  be  compared  across  different  scanners  and  provides  a  more  robust  value 
for  relative  comparison. 

In  summary,  in  this  proof  of  concept  study,  we  present  our  initial  clinical  experience  with 
prostate  RSI-MRI.  RSI-MRI  has  the  potential  to  substantially  improve  the  radiological  detection 
and  characterization  of  PCa.  Further  translational  studies  will  focus  on  defining  its  operating 
characteristics  with  respect  to  diagnosis  and  staging. 
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Figure  1.  MR  imaging  and  whole  mount  pathology  for  7  prostate  cancer  patients  with  EPE. 
Column  1 :  T2  -  weighted  images,  column  2:  Perfusion  Ktrans  maps,  column  3:  standard  apparent 
diffusion  coefficient  (ADC)  maps  (b  =  0,  100,  400,  and  800  s/mm2),  column  4:  RSI  z-score 
maps,  column  5:  whole  mount  pathology  with  tumor  and  area  of  EPE  identified.  Color  bar 
represents  zscores  from  0  to  7  for  the  RSI  maps. 


Figure  2.  MR  imaging  and  whole  mount  pathology  for  3  prostate  cancer  patients  without  EPE. 
Column  1:  T2  -  weighted  images,  column  2:  Perfusion  Ktrans  maps,  column  3:  standard  apparent 
diffusion  coefficient  (ADC)  maps  (b  =  0,  100,  400,  and  800  s/mm2),  column  4:  RSI  z-score 
maps,  column  5:  whole  mount  pathology  with  tumor  and  area  of  EPE  identified.  Color  bar 
represents  zscores  from  0  to  7  for  the  RSI  maps. 


Table  1.  MRI  scan  parameters  for  prostate  MRI  protocol  at  3T 

Pulse  sequence 

Parameters 

T2 

Axial  3D  TSE  T2  (Siemens  SPACE, 

TR/TE  3800-5040/101  ETL  13,  14  cm 

FOV,  256  x  256  matrix,  1.5  mm 
contiguous  slices) 

Diffusion- weighted  (Standard) 

(echoplanar,  TR/TE  3900/60,  21  x  26  cm 
FOV,  130  x  160  matrix,  3.6  mm  slices,  4 
NEX,  b-values  0,  100,  400,  800  s/nnn2) 

T 1  Dynamic  Perfusion  Imaging 

(Siemens  TWIST,  TR./TE  3. 9/ 1.4  ms,  12° 
flip  angle,  26  x  26  cm  FOV,  160  x  160 
matrix,  3.6  mm  slices,  4.75  s/acquisition 
over  6  minutes  with  15  s  injection  delay, 
image  analysis  using  iCAD  Versavue) 

Restriction  Spectrum  Imaging 

Spin  echo  EPI,  TR/TE  5500/137,  26  x  26 
cm  FOV,  128  x  96  matrix,  3.6  mm  slices, 

30  directions  at  each  b-value,  b-values  0, 

800,  1500,  4000  s/mm2. 

Table  2.  Radiologic  findings  -  documented  diagnostic  interpretations  Din  7  patients  with  biopsy-proven  localized 
prostate  cancer 

Case 

Location 

Axial  Clock 
Face  Location 

Capsular  Involvement 

Technical  Quality 

1 

Right  peripheral 
gland 

7-8  o’clock 

Focal  EPE  may  involve 
adjacent  right  neurovascular 
bundle.  Right  seminal 
vesicle  invasion. 

Excellent,  no  limitations 

2 

Right  central 
peripheral  basal 
mid  gland 

7  o’clock 

Intact,  no  suspicion  for 
involvement 

Excellent,  no  limitations 

3 

Right  peripheral 
basal  mid  gland 

8  o’clock 

Bulges  the  capsule  with  a 
broad  base  of  contact,  organ 
confined  disease 

Excellent,  no  limitations 

4 

Right  peripheral 
gland  lesion 

6-9  o’clock 

Bulges  the  capsule, 
suspicious  but  not  definitive 
for  extracapsular  extension 

Excellent,  no  limitations 

5 

R  peripheral 
gland  from  apex 
to  base 

5-10  o’  clock 

Blurring  and  irregularity  of 
right  capsule,  suspicious  for 
capsular  involvement 

Limited  spectral  quality,  does  not 
diminish  diagnostic  confidence 

6 

R  peripheral 
base  lesion 

7  o’  clock 

Bulges  the  capsule  with  a 
broad  base  of  attachment 

Excellent,  no  limitations 

7 

Right  peripheral 
mid  gland 

9-10  o’clock 

Bulges  the  capsule  laterally 

Flemorrhage  limits  sensitivity  and 
functional  characterization 

Table  3.  Pathologic  findings  determined  from  post  resection  whole  mount  pathology  in  7  patients 
undergoing  radical  prostatectomy 


Case 

Gleason 

Location  of  EPE 

Staging 

Age 

PSA 

Surgical  Margins 
at  Site  of  EPE 

1 

3+4 

Right  Posterior 

pT3bNxMx 

61 

5.8 

Negative 

2 

3+4 

Right  Posterior 

pT3aN0Mx 

55 

2.7 

Positive 

3 

4  +  5 

Right  Posterior 

pT3bpNlMx 

64 

8.2 

Negative 

4 

4  +  3 

Right  Posterior 

pT3aNlMx 

64 

4.2 

Positive 

5 

4  +  5 

Right  Posterior 

pT3aNlMx 

68 

7 

Positive 

6 

4  +  3 

Right  Posterior 

pT3apN0Mx 

56 

10 

Positive 

7 

4  +  3 

Right  Antero-lateral 

pT3apN0Mx 

67 

5.1 

Negative 

Table  4.  Quantitative  Diffusion  Imaging  -  range  of  ADC 
scores  and  corresponding  RSI  Z-scores 

Case 

ADC  (10A-3  mmA2/s) 

RSI  Z-score 

1 

1.08 

5.9 

2 

0.99 

5.7 

3 

1.01 

5.8 

4 

0.89 

4.6 

5 

1.26 

6.6 

6 

1.38 

4.2 

7 

1.12 

5.0 

MRI-derived  Restriction  Spectrum  Imaging  Cellularity  Index  is  Associated 
with  High  Grade  Prostate  Cancer  on  Radical  Prostatectomy  Specimens 
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Abstract 


Purpose:  We  evaluate  a  novel  magnetic  resonance  imaging  (MRI)  technique  to 
improve  detection  of  aggressive  prostate  cancer. 

Materials  and  Methods:  We  performed  a  retrospective  analysis  of  presurgical 
prostate  MRI  scans  using  an  advanced  diffusion  weighted  imaging  technique 
called  Restriction  Spectrum  Imaging  (RSI),  which  can  be  presented  as  a 
normalized  z-score  statistic.  Scans  were  acquired  prior  to  radical  prostatectomy. 
Prostatectomy  specimens  were  processed  using  whole  mount  sectioning  and 
regions  of  interest  (ROIs)  were  drawn  around  individual  prostate  cancer  (PCa) 
tumors.  Corresponding  ROIs  were  drawn  on  the  MRI  imaging  and  paired  with 
ROIs  in  regions  with  no  pathology.  RSI  z-score  and  conventional  apparent 
diffusion  coefficient  (ADC)  values  were  recorded  for  each  ROI.  Paired  t-test, 
ANOVA  and  logistic  regression  analyses  were  performed. 

Results:  We  evaluated  28  patients  with  64  regions  of  interest  (28  benign  and  36 
PCa).  The  mean  difference  in  RSI  z-score  (PCa  ROI  -  Benign  ROI)  was  2.17 
(SE  =  0.11;  p  <0.001 )  and  in  ADC  was  551  mm2/sec  (SE  =  80  mm2/sec;  paired  t- 
test,  p  <0.001).  The  differences  in  the  means  among  all  groups  (benign,  primary 
Gleason  3  and  primary  Gleason  4)  was  significant  for  both  RSI  z-score  (F3  64  = 
97.7,  p  <0.001 )  and  ADC  (F3,64  =  13.9,  p  <0.001 ).  A  t-test  was  performed  on  only 
PCa  tumor  ROIs  (n=36)  to  determine  prostate  cancer  aggressiveness  (Gleason  3 
vs.  Gleason  4)  revealing  that  RSI  z-score  was  still  significant  (p  =  0.03),  whereas, 
ADC  values  were  no  longer  significant  (p  =  0.08).  In  multivariable  analysis 
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adjusting  for  age  and  race,  RSI  z-score  was  associated  with  PCa  aggressiveness 
(OR  10.3,  95%CI:  1.4-78.0,  p=0. 02)  while  ADC  trended  to  significance  (p=0.07). 
Conclusions:  The  RSI  derived  normalized  cellularity  index  (RSI  z-score)  is 
associated  with  aggressive  prostate  cancer  as  determined  by  pathologic  Gleason 
scores.  Further  utilization  of  RSI  techniques  may  serve  to  enhance  standardized 
reporting  systems  for  PCa  in  the  future. 
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introduction 


One  current  focus  in  prostate  cancer  (PCa)  diagnosis  is  to  distinguish  indolent 
from  more  aggressive  disease  to  reduce  over-treatment.1  2  Magnetic  resonance 
imaging  (MRI)  may  be  a  non-invasive  imaging  biomarker  to  incorporate  into  PCa 
detection  and  treatment  strategies.3 

Recently,  MRI  has  been  increasingly  investigated  for  use  as  a  tool  in  the 
screening,  staging,  and  monitoring  of  prostate  cancer.4'6  Imaging  techniques 
such  as  diffusion-weighted  imaging  (DWI)  and  the  resultant  quantitative  apparent 
diffusion  coefficient  (ADC)  have  shown  correlation  with  PCa;  however,  they  have 
significant  limitations  regarding  tumor  conspicuity  and  localization.7  8  A  novel, 
advanced  diffusion-based  imaging  technique,  called  Restriction  Spectrum 
Imaging  (RSI),  has  been  modified  from  previous  studies  in  brain  cancer  detection 
to  apply  to  patients  with  prostate  cancer.8, 9  Some  benefits  of  the  RSI  technique 
is  to  provide  reduced  spatial  distortion,  enhanced  tumor  contrast-to-noise  over 
conventional  diffusion  weighted  imaging  (DWI)  and  a  normalized  in  vivo  measure 
of  cellularity. 

Within  individual  tumor  regions  of  interest,  we  investigate  the  association  of  final 
pathologic  Gleason  score  from  whole-mount  prostatectomy  specimens  with  the 
RSI  cellularity  index  as  compared  to  the  current  standard,  ADC.10'16  Our  primary 
outcome  is  the  detection  of  primary  pattern  Gleason  >4  prostate  cancer. 
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Materials  and  Methods 


Patients 

All  patients  were  previously  diagnosed  with  prostate  cancer  via  standard 
transrectal  ultrasound  guided  prostate  biopsy  after  prostate  specific  antigen 
(PSA)  elevation  or  abnormal  digital  rectal  examination  (DRE).  Pre-surgical  pelvic 
MRI  is  routinely  used  at  our  institution  to  identify  extraprostatic  extension  (EPE) 
in  order  to  guide  nerve-sparing  surgery.  Data  were  collected  from  chart  review. 

MRI  and  RSI 

Patients  underwent  standard  T2,  perfusion  (with  Gadolinium)  and  diffusion 
protocols  at  3T  (Siemens,  Erlangen  Germany)  with  an  endorectal  coil  prior  to 
radical  prostatectomy.  Table  1  shows  pulse  sequence  details.  The  ADC  maps 
used  to  draw  ROIs  were  generated  from  the  low  b-value,  800s/mm2,  derived 
from  the  same  spectrum  of  b-values  used  in  the  RSI  protocol.  ADC  maps  were 
corrected  for  spatial  distortion  (ref  17).  The  restriction  spectrum  diffusion  tensor 
imaging  protocol  parameters  include  b-values  of  0,  800,  1500,  4000  s/mm2  in  30 
unique  diffusion  directions  for  each  nonzero  b-value.  Imaging  time  for  this 
sequence  is  approximately  3.5  minutes.  RSI  cellularity  maps  were  reconstructed 
based  on  all  b-values8.  The  RSI  cellularity  maps  were  then  standardized  across 
the  sample,  using  the  mean  and  standard  deviations  of  normal  prostate  signal 
from  the  raw  RSI  maps  in  20  patients  to  produce  z-score  maps.  RSI  maps  were 
also  corrected  for  spatial  distortion.17 
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Pathology 

After  prostatectomy,  whole  mount  histopathology  was  performed  on  4  pm  thick 
sections  of  each  specimen.  A  Gleason  score  was  assigned  to  each 
representative  tumor  location.  If  two  tumors  were  located,  the  Gleason  score  for 
each  was  assessed  independently.  The  histopathology  was  evaluated  and  the 
boundaries  of  tumor  versus  benign  tissue  were  identified  by  a  uropathologist. 

Outcomes 

We  defined  our  primary  outcome  as  pathologic  primary  Gleason  score  of  4, 
which  means  that  Gleason  4  is  the  dominant  histologic  architecture  and  includes 
4+3,  4+4,  and  4+5  Gleason  patterns.  The  pathologic  Gleason  score  is  currently 
the  standard  of  reference  for  PCa  aggressiveness.  Additionally,  the  ability  of 
imaging  to  detect  secondary  Gleason  patterns  may  be  minimal;  therefore,  herein 
we  focus  on  primary  Gleason  patterns.  Our  primary  predictor  variable  was  the 
normalized  cellularity  index  called  the  “RSI  z-score”.  The  most  commonly  utilized 
tool  to  identify  and  classify  aggressive  cancer  on  MRI  currently  is  the  ADC  value 
from  DWI;  therefore,  the  RSI  z-score  was  compared  with  ADC  to  assess  the 
predictive  value  in  differentiating  cancer  from  normal  regions  of  interest  (ROI). 

Each  patient  had  at  least  one  identified  region  of  cancerous  tissue.  If  two  areas 
of  cancer  were  detected,  each  region  was  evaluated  and  assigned  a  separate 
Gleason  score,  ADC,  and  z-score.  Tumor  ROI’s  were  drawn  based  on  the 
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pathology  in  combination  with  ADC  images  that  had  been  corrected  for  spatial 
distortion.  A  benign  region  of  interest  (ROI)  was  defined  in  a  region  of  the 
prostate  seen  to  be  free  of  prostate  cancer  on  the  whole  mount  histology.  ADC 
and  RSI  z-score  values  were  recorded  for  all  ROIs. 

Statistical  Analysis 

Each  patient  had  at  least  one  ROI  of  cancer  and  one  ROI  of  benign  tissue. 
Correlation  between  the  RSI-derived  z-score  and  ADC  was  determined  by  a 
Pearson  correlation  test  assuming  normal  distribution.  In  order  to  investigate  the 
association  of  RSI  z-score  and  primary  Gleason  pattern  4  PCa  vs.  pattern  3  PCa, 
a  t-test  test  was  performed.  In  order  to  compare  the  utility  of  MRI  techniques  (RSI 
z-score  vs.  ADC)  for  detecting  aggressive  cancer,  we  compared  ROIs 
representing  pathologically  benign  tissue  with  those  representing  increasing 
aggressive  prostate  cancer  (benign  vs.  Gleason  3  vs.  Gleason  4  primary 
patterns)  using  ANOVA  analysis  (F-test).  After  removing  the  values  for  the 
benign  ROIs,  we  also  assessed  variation  in  MRI  values  among  different  grades 
of  cancerous  tissue  aggressiveness  by  performing  a  comparative  t-test. 
Multivariable  analysis  included  an  ordinal  logistic  regression  (benign  vs.  Gleason 
3  vs.  Gleason  4)  and  binary  logistic  regression  (Gleason  3  vs.  Gleason  4).  p- 
values  less  than  0.05  were  considered  statistically  significant  using  the  statistical 
package  SPSS  v.21  (IBM,  Chicago,  IL.). 
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Results 


We  identified  28  patients  who  underwent  preoperative  MRI  with  RSI  and 
subsequent  whole  mount  pathology  after  radical  prostatectomy,  with  surgery 
taking  place  between  May  2012  and  May  2013.  Demographics  are  displayed  in 
Table  2.  Figure  1  shows  representative  examples  of  RSI  across  the  Gleason 
spectrum  and  illustrates  data  collection  methods.  White  arrows  within  the  higher 
grade  RSI  maps  show  areas  of  signal  void  that  could  be  interpreted  as  false 
positives  on  the  corresponding  ADC  maps.  We  identified  64  regions  of  interest 
(28  benign  and  36  PCa).  Eight  patients  had  two  distinct  tumor  ROIs  within  one 
specimen.  Seven  of  those  patients  had  discordant  tumors  (two  different  Gleason 
scores)  and  only  one  patient  had  two  concordant  tumor  ROIs  (both  Gleason 
3+4).  The  RSI  z-score  data  for  all  ROIs,  grouped  by  Gleason  score,  is  shown  in 
Figure  2. 

ADC  and  RSI  z-scores  are  normally  distributed  and  did  not  need  further 
transformation.  There  is  a  significant  correlation  between  ADC  and  RSI  z-score 
(Pearson  R  =  0.69;  n=64,  p  <0.001 )  (see  supplementary  figure).  ADC  and  RSI  z- 
score  are  collinear  (collinearity  index  8.4);  therefore,  they  cannot  be  placed  within 
the  same  multivariate  analytic  model. 

The  mean  cellularity  index  (RSI  z-score)  in  PCa  was  2.53  (SE  =  0.10)  and  in 
benign  tissue  0.39  (SE  =  0.12;  p<0.001 ).  The  mean  ADC  for  PCa  was  1169 
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mm2/sec  (SE  =  67  mm2/sec)  and  benign  was  1679  mm2/sec  (SE  =  81  mm2/sec,  p 
<0.001). 

For  the  univariate  and  multivariate  analyses,  the  data  were  grouped  by  primary 
Gleason  pattern,  either  primary  Gleason  3  (3+3  or  3+4)  or  Gleason  4  (4+3,  4+4, 
or  4+5).  In  univariable  analysis,  we  investigate  the  association  of  RSI  z-score  and 
ADC  with  PCa  on  pathology  (benign  vs.  Gleason  3  vs.  Gleason  4)  using  ANOVA 
analysis.  Both  were  able  to  distinguish  benign  from  increasingly  malignant 
prostate  cancer  (both  p<0.001)  (Table  3).  Figure  3  displays  a  box-plot  developed 
from  individual  ROI  RSI  z-scores  corresponding  to  benign  tissue  primary 
Gleason  3  pattern  or  primary  Gleason  4  pattern  PCa  tumors.  When  removing  the 
benign  regions  of  interest  as  to  only  compare  low-grade  (primary  Gleason  3)  to 
high-grade  PCA  (primary  Gleason  4),  RSI  technique  was  able  to  distinguish  the 
two  groups  (p=0.03)  and  ADC  trended  to  significance  (p=0.08).  We  then 
investigated  these  same  parameters  in  multivariable  analysis  and  determined 
that  both  RSI  z-score  and  ADC  were  able  to  distinguish  between  the  three 
groups  in  ordinal  regression  analysis  adjusting  for  age  and  race  (white  vs.  non¬ 
white);  though  the  log  odds  of  the  parameter  estimates  suggest  an  improved 
distinction  of  the  groups  by  RSI  (RSI  z-score  and  ADC,  p<0.001).  However,  when 
determining  the  difference  in  the  detection  of  low-grade  and  high-grade  PCa,  a 
higher  RSI  z-score  was  significantly  associated  with  the  higher-grade  primary 
Gleason  4  pattern  (Odds  ratio  10.3  (1.4-78.0;  p=0.02)  and  ADC  showed  a  trend 
in  distinguishing  between  the  two  Gleason  patterns  (p=0.07)  (Table  4). 
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Discussion 


Magnetic  Resonance  Restriction  Spectrum  Imaging  normalized  cellularity  index 
(RSI  z-score)  is  able  to  distinguish  aggressive  PCa  (primary  Gleason  score  of  4 
compared  to  3)  in  our  population  of  men  undergoing  radical  prostatectomy. 
Importantly,  the  RSI  technique  has  displayed  at  least  similar  ability  to  distinguish 
Gleason  grade  to  the  current  reference  standard,  ADC  values. 

Multiple  studies  have  described  the  ability  of  ADC  to  detect  PCa.  However,  the 
distinction  between  prostate  cancer  aggressiveness  has  been  less  investigated. 
Donati  et  al.  found  that  mean  ADC  could  distinguish  Gleason  6  from  7+  tumors  in 
131  men  undergoing  prostatectomy  (AUC  0.706)  and  in  another  paper  discussed 
the  use  of  10th  percentile  ADC  correlation  to  aggressiveness.10, 16  The  10th 
percentile  ACD  was  also  used  in  combination  with  mean  ADC,  T2-weighted 
skewness  and  k  trans  to  distinguish  PCa  using  computer  aided  diagnosis.18 
Moreover,  a  recent  study  has  suggested  that  ADC  entropy  rather  than  mean 
ADC  could  better  discriminate  the  proportion  of  Gleason  4  cancer  among 
Gleason  3+4  and  4+3  tumors.19  The  distinction  in  primary  Gleason  pattern  may 
have  significant  clinical  implications  regarding  prostate  cancer  management 
decision-making. 
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Primary  Gleason  4  pattern  is  an  aggressive  cancer  with  patients  experiencing 
higher  rates  of  biochemical  (PSA)  failure  after  prostatectomy,  systemic 
recurrence  and  PCa  mortality.20  Therefore,  knowledge  of  high-grade  cancer  prior 
to  making  management  decisions  would  be  in  helpful  determining  treatment 
strategy.  For  example,  men  without  Gleason  4  pattern  PCa  are  more  ideal 
candidates  for  active  surveillance.22  Serial  imaging  may  indicate  progression  of 
disease,  assisting  urologists  in  deciding  when  a  biopsy  off  of  protocol  is 
warranted.  Additionally,  PCa  grade  may  influence  the  urologic  surgeon  to 
perform  a  pelvic  lymph  node  dissection  at  the  time  of  prostatectomy  due  to 
increased  risk  of  nodal  disease.21 

Currently,  ADC  serves  as  the  most  discriminatory  parameter  to  assist  radiologists 
for  the  detection  of  cancer.  Moreover,  recent  examination  of  ADC  and  prostate 
cancer  has  shown  the  association  of  ADC  and  PCa  aggressiveness.10, 16  While 
we  do  show  that  ADC  can  differentiate  the  presence  of  cancer  or  not,  our  study 
shows  ADC  is  less  able  to  determine  the  subtlety  of  primary  pattern  Gleason  3 
vs.  Gleason  4  PCa.  Possible  reasons  include  the  proportion  of  pattern  3  vs.  4 
disease  in  our  population  compared  to  prior  populations,  the  b-values  used  in 
determining  the  ADC,  how  the  ROI  was  chosen,  the  amount  of  stromal  reaction, 
and  technical  factors  such  as  degree  of  hemorrhage.  The  differences  between 
ADC  and  RSI  z-score  in  our  study  are  small  however. 
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RSI  techniques  offer  advantages  when  compared  to  conventional  DWI  and  ADC 
maps.  For  example,  one  of  the  challenges  of  standard  diffusion  imaging,  is  that 
the  ADC  values  are  not  standardized  across  MRI  scanners.  The  z-score  is  a 
standardized  statistical  method  and  inherently  normalizes  across  the  patient 
pool.  Thus,  the  RSI  z-score  is  a  value  that  could  potentially  be  compared  across 
different  scanners  and  institutions  and  provide  a  more  robust  value  for  relative 
comparison. 

In  addition,  conventional  DWI/ADC  suffers  from  geometric  distortion  and  can  be 
difficult  to  interpret  by  clinicians  and  untrained  radiologists.  Distortion  correction 
techniques  previously  optimized  in  the  brain  for  GBM,  are  incorporated  into  our 
restriction  spectrum  imaging  post-processing  stream  in  order  to  derive  spatially 
corrected  cellularity  maps.  Distortion  correction  techniques  are  not  routinely 
employed  in  conventional  DWI  and  resultant  ADC  maps.  Thus,  the  RSI  maps 
can  be  co-registered  with  T2  weighted  anatomic  images  with  voxel  accuracy. 

This  has  potential  implications  for  more  accurate  detection  of  extraprostatic 
extension  (ref  Rebecca’s  accepted  paper)  and  more  accurate  MRI-fused 
ultrasound  targeted  biopsy  results. 

Because  of  its  greater  sensitivity  to  restriction  rather  than  hindered  diffusion,  RSI 
may  be  less  subject  to  hemorrhage,  inflammatory  processes,  and  benign  nodules 
in  the  transitional  zone,  all  of  which  can  exhibit  lower  ADC  values  leading  to  false 
positives.  Theoretically,  RSI-MRI  reduces  extracellular  signal  by  focusing  on  the 
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signal  emanating  from  intracellular  tumor  cells  (restriction  diffusion)  and  less  from 
the  extracellular  signal  (hindered  diffusion). 9  This  will  need  to  be  rigorously 
tested  in  future  ROC  performance  studies. 

In  addition,  in  regions  where  there  is  overt  chemical  dephasing  from  gross 
calcium  or  hemorrhage  or  other  etiology,  resulting  in  signal  void,  ADC  maps  will 
show  these  areas  as  low  signal,  which  can  be  interpreted  inaccurately,  leading  to 
a  false  positive.  White  arrow  heads  in  Figure  1  show  two  such  examples.  For 
example,  in  the  Gleason  5+4  case,  the  anterior  region  of  signal  void  shows  up  as 
dark  (low)  on  the  ADC  maps,  potentially  a  false  positive,  while  in  the  RSI  maps, 
this  is  clearly  interpreted  as  an  area  of  signal  void,  not  tumor.  Thus,  RSI  offers  a 
number  of  potential  advantages  when  compared  to  DWI/ADC. 

Certain  limitations  of  this  study  include  small  sample  size  and  retrospective  data 
collection.  We  have  overcome  the  small  sample  size  by  using  each  patient  as  his 
own  control  to  provide  a  paired  analysis  by  using  known  benign  tissue.  However, 
because  of  the  small  sample  size  we  have  fewer  patients  with  extremely  high- 
grade  cancer  (Gleason  5)  and  low-grade  (Gleason  6  or  less);  therefore,  we 
dichotomized  based  on  the  primary  Gleason  pattern.  Additionally,  we  feel 
imaging  is  unlikely  to  visualized  smaller  amounts  (secondary  patterns)  of  PCa 
architecture.  We  do  have  a  selection  bias  regarding  our  patient  population  as  all 
patients  underwent  radical  prostatectomy.  Therefore,  our  results  may  not 
necessarily  apply  to  patients  in  the  general  PSA  screening  population  undergoing 
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prostate  biopsy  is  an  area  of  future  investigation.  The  lack  of  patients  in  this 
study  with  pathologies  at  the  extremes  shows  the  need  for  a  broader  study. 

Conclusions 

Magnetic  Resonance  Restriction  Spectrum  Imaging,  which  can  provide  a 
normalized  index  of  cellularity  (RSI  z-score),  is  associated  with  detection  of 
aggressive  PCa  as  defined  by  Gleason  score.  Additionally,  RSI  includes 
correction  of  spatial  distortion,  a  normalized  measure  of  cellularity  and  in  general 
increased  conspicuity  when  compared  to  conventional  DWI/ADC.  RSI 
technology  warrants  prospective  evaluation  in  the  PCa  diagnostic  arena. 


15 


References 

1 .  Johansson,  J.  E.,  Andren,  O.,  Andersson,  S.  O.  et  al.:  Natural  history  of 
early,  localized  prostate  cancer.  JAMA,  291:  2713,  2004 

2.  Loeb,  S.,  Bjurlin,  M.  A.,  Nicholson,  J.  et  al.:  Overdiagnosis  and 
overtreatment  of  prostate  cancer.  Eur  Urol,  65:  1046,  2014 

3.  van  den  Bergh,  R.  C.,  Ahmed,  H.  U.,  Bangma,  C.  H.  et  al.:  Novel  tools  to 
improve  patient  selection  and  monitoring  on  active  surveillance  for  low-risk 
prostate  cancer:  a  systematic  review.  Eur  Urol,  65:  1023,  2014 

4.  Chamie,  K.,  Sonn,  G.  A.,  Finley,  D.  S.  et  al.:  The  role  of  magnetic 
resonance  imaging  in  delineating  clinically  significant  prostate  cancer. 
Urology,  83:  369,  2014 

5.  Park,  B.  H.,  Jeon,  H.  G.,  Jeong,  B.  C.  et  al.:  Influence  of  Magnetic 
Resonance  Imaging  in  the  Decision  to  Preserve  or  Resect  Neurovascular 
Bundles  at  Robotic  Assisted  Laparoscopic  Radical  Prostatectomy.  J  Urol, 
2014 

6.  Lee,  D.  J.,  Ahmed,  H.  U.,  Moore,  C.  M.  et  al.:  Multiparametric  magnetic 
resonance  imaging  in  the  management  and  diagnosis  of  prostate  cancer: 
current  applications  and  strategies.  Curr  Urol  Rep,  15:  390,  2014 

7.  Hambrock,  T.,  Somford,  D.  M.,  Huisman,  H.  J.  et  al.:  Relationship  between 
apparent  diffusion  coefficients  at  3.0-T  MR  imaging  and  Gleason  grade  in 
peripheral  zone  prostate  cancer.  Radiology,  259:  453,  201 1 

8.  White,  N.  S.,  McDonald,  C.  R.,  Farid,  N.  et  al.:  Improved  conspicuity  and 
delineation  of  high-grade  primary  and  metastatic  brain  tumors  using 
"restriction  spectrum  imaging":  quantitative  comparison  with  high  B-value 
DWI  and  ADC.  AJNR  Am  J  Neuroradiol,  34:  958,  2013 

9.  White,  N.  S.,  Leergaard,  T.  B.,  D'Arceuil,  H.  et  al.:  Probing  tissue 
microstructure  with  restriction  spectrum  imaging:  Histological  and 
theoretical  validation.  Hum  Brain  Mapp,  34:  327,  2013 

10.  Donati,  O.  F.,  Afaq,  A.,  Vargas,  H.  A.  et  al.:  Prostate  MRI:  evaluating 
tumor  volume  and  apparent  diffusion  coefficient  as  surrogate  biomarkers 
for  predicting  tumor  Gleason  score.  Clin  Cancer  Res,  20:  3705,  2014 

1 1 .  Kobus,  T.,  Vos,  P.  C.,  Hambrock,  T.  et  al.:  Prostate  cancer 
aggressiveness:  in  vivo  assessment  of  MR  spectroscopy  and  diffusion- 
weighted  imaging  at  3  T.  Radiology,  265:  457,  2012 

12.  Oto,  A.,  Yang,  C.,  Kayhan,  A.  et  al.:  Diffusion-weighted  and  dynamic 
contrast-enhanced  MRI  of  prostate  cancer:  correlation  of  quantitative  MR 


16 


parameters  with  Gleason  score  and  tumor  angiogenesis.  AJR  Am  J 
Roentgenol,  197:  1382,  2011 

13.  Turkbey,  B.,  Shah,  V.  P.,  Pang,  Y.  et  al.:  Is  apparent  diffusion  coefficient 
associated  with  clinical  risk  scores  for  prostate  cancers  that  are  visible  on 
3-T  MR  images?  Radiology,  258:  488,  201 1 

14.  Vargas,  H.  A.,  Akin,  O.,  Franiel,  T.  et  al.:  Diffusion-weighted  endorectal 
MR  imaging  at  3  T  for  prostate  cancer:  tumor  detection  and  assessment  of 
aggressiveness.  Radiology,  259:  775,  2011 

15.  Vargas,  H.  A.,  Donati,  O.  F.,  Wibmer,  A.  et  al.:  Association  Between 
Penile  Dynamic  Contrast-Enhanced  MRI-Derived  Quantitative  Parameters 
and  Self-Reported  Sexual  Function  in  Patients  with  Newly  Diagnosed 
Prostate  Cancer.  J  Sex  Med,  2014 

16.  Donati,  O.  F.,  Mazaheri,  Y.,  Afaq,  A.  et  al.:  Prostate  cancer 
aggressiveness:  assessment  with  whole-lesion  histogram  analysis  of  the 
apparent  diffusion  coefficient.  Radiology,  271:  143,  2014 

17.  Holland,  D.,  Kuperman,  J.  M.,  Dale,  A.  M.:  Efficient  correction  of 
inhomogeneous  static  magnetic  field-induced  distortion  in  Echo  Planar 
Imaging.  Neuroimage,  50:  175,  2010 

18.  Peng,  Y.,  Jiang,  Y.,  Antic,  T.  et  al.:  Validation  of  quantitative  analysis  of 
multiparametric  prostate  MR  images  for  prostate  cancer  detection  and 
aggressiveness  assessment:  a  cross-imager  study.  Radiology,  271:  461, 
2014 

19.  Rosenkrantz,  A.  B.,  Triolo,  M.  J.,  Melamed,  J.  et  al.:  Whole-lesion 
apparent  diffusion  coefficient  metrics  as  a  marker  of  percentage  Gleason 
4  component  within  Gleason  7  prostate  cancer  at  radical  prostatectomy.  J 
Magn  Reson  Imaging,  2014 

20.  Tollefson,  M.  K.,  Leibovich,  B.  C.,  Slezak,  J.  M.  etal.:  Long-term 
prognostic  significance  of  primary  Gleason  pattern  in  patients  with 
Gleason  score  7  prostate  cancer:  impact  on  prostate  cancer  specific 
survival.  J  Urol,  175:  547,  2006 

21 .  Allaf,  M.  E.,  Partin,  A.  W.,  Carter,  H.  B.:  The  importance  of  pelvic  lymph 
node  dissection  in  men  with  clinically  localized  prostate  cancer.  Rev  Urol, 
8:112,2006 

22.  Klotz,  L.,  Zhang,  L.,  Lam,  A.  et  al.:  Clinical  results  of  long-term  follow-up  of 
a  large,  active  surveillance  cohort  with  localized  prostate  cancer.  J  Clin 
Oncol,  28:  126,  2010 


17 


18 


Tables 


Table  1 :  MRI  scan  parameters  for  prostate  MRI  protocol  at  3T. 

Table  2:  Demographics:  Perioperative  demographics  for  28  patients  who 
underwent  MRI  with  an  endorectal  coil  and  subsequently  underwent  radical 
prostatectomy. 

Table  3:  Univariable  analysis:  Restriction  spectrum  imaging  (RSI)  z-score  and 
apparent  diffusion  coefficient  (ADC)  are  independently  investigated  comparing 
individual  regions  of  interest  to  their  corresponding  primary  Gleason  pattern.  The 
top  half  of  the  table  uses  ANVOA  analysis  to  determine  differences  in  detection 
comparing  benign  vs.  primary  Gleason  3  vs.  Primary  Gleason  4  prostate  cancer. 
The  bottom  half  of  the  table  only  compares  Gleason  3  vs.  4  cancers  with  a  t-test. 

Table  4:  Multivariable  analysis:  Restriction  spectrum  imaging  (RSI)  z-score  and 
apparent  diffusion  coefficient  (ADC)  are  independently  investigated  comparing 
individual  regions  of  interest  to  their  corresponding  primary  Gleason  pattern.  The 
top  half  of  the  table  uses  ordinal  logistic  regression  analysis  to  determine 
differences  in  detection  comparing  benign  vs.  primary  Gleason  3  vs.  Primary 
Gleason  4  prostate  cancer  (Log  Odds).  The  bottom  half  of  the  table  only 
compares  Gleason  3  vs.  4  cancers  to  determine  with  binary  logistic  regression 
(Odds  Ratio). 
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Figures 


Figure  1 :  Representative  images  showing  RSI  z-score  maps  across  Gleason 
scores:  The  y-axis  shows  the  pathologic  Gleason  score  with  the  x-axis 
designating  the  MRI  sequence.  The  last  column  displays  the  whole  mount 
pathology  with  the  corresponding  cancer  region  of  interest  circled  in  black.  The 
star  in  the  top  right  pathologic  figure  represents  the  pattern  3+3  prostate  cancer 
while  the  other  lesions  are  4+3.  White  arrowheads  in  the  higher  grade  patients 
show  areas  of  signal  void  which  could  be  interpreted  as  false  positives  on  the 
ADC  maps. 

Figure  2:  RSI-z  Score  Value  Grouped  by  Pathologic  Gleason  Score:  The  y  -axis 
represents  the  RSI  z-score  derived  from  a  given  region  of  interest.  The  x-axis 
demonstrates  the  pathological  Gleason  scores  in  increasing  levels  of 
aggressiveness  from  left  to  right.  Each  data  point  represents  one  region  of 
interest  corresponding  to  a  location  on  the  whole  mount  prostatectomy  specimen 
contoured  by  a  GU  pathologist. 

Figure  3:  Box  Plot  of  RSI  z-Score  for  Primary  Gleason  Pattern:  The  box-plot 
represents  the  RSI  z-Score  for  benign,  pathologic  primary  pattern  Gleason  3,  or 
pathologic  primary  pattern  Gleason  4  prostate  cancer. 
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Table  1 


MRI  scan  parameters  for  prostate  MRI  protocol  at  3T 

Pulse  sequence 

Parameters 

T2 

Axial  3D  TSE  T2  (Siemens  SPACE, 

TR/TE  3800-5040/101  ETL  13,  14  cm 

FOV,  256  x  256  matrix,  1.5  mm 
contiguous  slices) 

Diffusion-weighted  (Standard) 

(echoplanar,  TR/TE  3900/60,  21  x  26  cm 
FOV,  130  x  160  matrix,  3.6  mm  slices,  4 
NEX,  b-values  0,  100,  400,  800  s/mm2) 

T 1  Dynamic  Perfusion  Imaging 

(Siemens  TWIST,  TR./TE  3. 9/ 1.4  ms,  12° 
flip  angle,  26  x  26  cm  FOV,  160  x  160 
matrix,  3.6  mm  slices,  4.75  s/acquisition 
over  6  minutes  with  15  s  injection  delay, 
image  analysis  using  iCAD  Versavue) 

Restriction  Spectrum  Imaging 

Spin  echo  EPI,  TR/TE  5500/137,  26  x  26 
cm  FOV,  128  x  96  matrix,  3.6  mm  slices, 

30  directions  at  each  b-value,  b-values  0, 

800,  1500,  4000  s/mm2. 
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Table  2 


ID  Age  Race 

Body 

Mass 

Index 

Preoperative 

PSA 

Biopsy 

Gleason 

Prostatectomy 

Gleason 

Pathologic 

Stage 

PIRADS 

Score 

1  63 

White 

20.1 

7 

3+4 

3+3 

PT2 

4 

2  53 

Asian 

16.9 

5.5 

3+4 

3+4 

pT2  Nx  Mx 

4 

3  58 

White 

32.7 

3.8 

3+4 

3+4 

pT2c  Nx  Mx 

3 

4  44 

White 

26.6 

2.8 

3+4 

3+4 

pT2a  NO  Mx 

3 

5  53 

White 

22.7 

3.9 

4+3 

3+4 

pT3a  pNO 

Mx 

2 

6  50 

White 

31.2 

8.9 

3+4 

3+4 

pT2  Nx  Mx 

4 

7  65 

Unknown  29 

5.4 

3+4 

3+4 

pT2  Nx  Mx 

2 

8  61 

White 

29.8 

6.6 

4+3 

3+4 

pT2  NO  Mx 

4 

9  60 

- 

25.5 

6 

3+3 

3+4 

pT2  pNx  Mx 

3 

1061 

Other 

26 

3.4 

3+4 

3+4 

pT2c  NO  Mx 

3 

11  64 

White 

23.1 

5.8 

3+4 

3+4 

pT2c  Nx  Mx 

2 

1255 

White 

27.5 

4.7 

3+4 

3+4 

pT2b  NO  Mx 

3 

13  68 

White 

- 

6.7 

3+4 

3+4 

pT2c  NO  Mx 

3 

1461 

White 

21.4 

9.2 

3+3 

3+4 

pT2  NO  Mx 

4 

15  59 

White 

31.4 

4.43 

3+3 

3+4 

pT2c  Nx  Mx 

3 

1671 

White 

29.7 

7.3 

3+3 

3+4 

PT2 

3 

17  55 

White 

29.3 

2.8 

3+5 

3+4 

pT3a 

3 

1861 

White 

31.4 

5.8 

4+3 

3+4 

pT3b 

5 

1971 

White 

28.4 

12.8 

4+3 

4+3 

pT3a  pNI 

Mx 

5 

20  67 

White 

30.4 

5.1 

3+4 

4+3 

pT3a  pNO 

Mx 

4 

21  56 

Other 

28 

- 

3+4 

4+3 

pT3a  NO  Mx 

4 

22  64 

White 

26.9 

4.2 

3+5 

4+3 

pT3a 

5 

23  65 

White 

38.5 

5.8 

4+3 

4+3 

pT3a 

5 

24  62 

Black 

25.1 

4.8 

4+3 

4+3 

pT2a  NO  Mx 

0 

25  45 

Black 

23.5 

6.5 

4+3 

4+4 

pT2c  NO  Mx 

1 

26  70 

Filipino 

35.5 

7.5 

3+4 

4+5 

pT3a  pNx 

Mx 

3 

27  62 

White 

25.8 

4.6 

4+5 

4+5 

pT2  NO  Mx 

4 

28  65 

White 

30.1 

8.2 

5+4 

4+5 

pT3b 

4 
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Table  3 


Sample 

Statistical 

P 

Size 

Method 

Means 

Test 

Value 

Detection  of  Increasingly  Aggressive 

Mean  (Standard 

Cancer* 

Error) 

F  test 

Restricted  Spectrum  Imaging 

64 

(RSI)  z-score 

Benign 

0.65  (0.12) 

97.7 

<0.001 

Primary 
Gleason  3 

2.4  (0.61) 

Primary 
Gleason  4 

2.9  (0.51) 

Apparent  Diffusion 

64 

Coefficient  (ADC) 

Benign 

1680  (428) 

13.9 

<0.001 

Primary 
Gleason  3 

1237  (425) 

Primary 
Gleason  4 

967  (221) 

Prostate  Cancer  Aggressiveness** 

Mean  Difference 

t  test 

Restricted  Spectrum  Imaging 

36 

(RSI)  z-score 

Apparent  Diffusion 

34 

0.5 

2.22 

0.033 

36 

Coefficient  (ADC) 

34 

269 

1.81 

0.079 
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Table  4 


Parameter 

Coefficient3 

Sample 

(95%  Confidence 

P 

McFadden's 

Size 

Method 

Interval) 

Value 

D 

Detection  of  Increasingly 

Aggressive  Prostate  Cancel 

Restricted  Spectrum 

64 

Imaging  (RSI)  z-score 

Benign 

0 

<0.001 

0.613 

Primary 
Gleason  3 

9.5  (-1.9-17.0) 

0.014 

Primary 
Gleason  4 

15.8  (6.4-25.2) 

0.001 

Apparent  Diffusion 

64 

Coefficient  (ADC) 

Benign 

0 

<0.001 

0.202 

Primary 
Gleason  3 

1.7  (-3.2  -  6.7) 

0.483 

Primary 
Gleason  4 

4.5  (-0.5  -  9.6) 

0.076 

Prostate  Cancer  Aggressiveness0 

Wald  test 

Odds  Ratio 
(95%  Confidence 
Interval) 

P 

Value 

36 

Restricted  Spectrum 
Imaging  (RSI)  z-score 

5.1 

10.3  (1.4-78.0) 

0.024 

36 

Apparent  Diffusion 
Coefficient  (ADC) 

3.3 

0.9  (0.9-1.00) 

0.069 
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Figure  1 
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Figure  3 


Restricted  Spectrum  Imaging  MRI  Z-Score  (Cellularity  Index) 


